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Abstract 

This paper deals with two separate but related results. First we consider weak 
solutions to a parabolic operator with Hormander vector fields. Adapting the iter- 
ation scheme of Jurgen Moser for elliptic and parabolic equations in W n we show 
a parabolic Harnack inequality. Then, after proving the Harnack inequality for 
weak solutions to equations of the form u t = Y2^i( a ij^j u ) we use this to show 
Holder continuity. We assume the coefficients are bounded and elliptic. The itera- 
tion scheme is a tool that may be adapted to many settings and we extend this to 
nonlinear parabolic equations of the form u t = —X*Aj(Xju). With this we show 
both a Harnack inequality and Holder continuity of weak solutions. 
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Notation Use 



the n-dimensional Euclidean space 

the product of the n-dimensional Euclidean space and a time 
interval (0,T) 

x a point in n-dimensional Euclidean space, x = {xi,x 2 , ■■■X n } T 

d cc (x,y) the Carnot-Caratheodory (CC) distance from x to y. (see pg. 
20) 

B cc (x,r) a CC ball centered at x with radius r. B cc (x,r) = {y e 
R n | d cc (x,y) < r}. Also B(x,r), B r . (see pg. 20) 

\B cc (x, r)\ the volume of the ball. \B cc (x, r) \ — J B dx. Here dx is Lebesgue 
measure 

R an open bounded rectangle in R™ +1 . It is the product of a 

Carnot-Caratheodory ball and a time interval: -B cc (x, r) x (0, T) 

u(x,t) a solution to a partial differential equation. u(x,t) : IR n+1 — >■ M 

Xi a vector field. 

a collection of vector fields. 

XiU the derivative of w in the direction of JQ. 

X* the formal adjoint of X;. (see pg. 18) 

Vw the spatial gradient of u. 

|Vn| the norm of the spatial gradient of n. |Vn| = 

Er=i(^/^) 2 ] i/2 

Von the horizontal gradient of n. Vo - " = YlT=i XiuXi. (see pg. 19) 

|Vo«| the norm of the horizontal gradient of u. |Von| = 

[E™i(^)T /2 

{aij(x, t)} an m x m matrix of whose coefficients are functions of x and t 
L 2 (R) the Lebesgue space of functions on R. It consists of functions 

u for which J R u 2 dxdt < oo. (see pg. 23) 
L°°(R) consists of all functions w for which ess sup w < oo 



Notation Use 

W 1,2 (R) the Sobolev space of functions on R. It consists of functions 
u E L 2 (R) with the weak derivatives of u in L 2 (R). 

S 1,2 (R) the Sub-elliptic Sobolev space of functions on R. It consists of 
functions u G L 2 (R) with the weak derivatives of u in L 2 (R). 
(see pg. 23) 

u t the derivative of u in time. u t = du/dt 



0.1 Introduction 



l 



This paper deals with two versions of the heat equation. The classical heat equation 
tells us the temperature at a given location in a material at a given time. It is simply 
modeled by 

u t = Au 

in some region. The derivatives here are the usual ones taken along the coordinates 
axes. This equation has been studied by many people and is very well known. This 
equation may also be written as: 

u t = div (Vtt). 

A typical property of interest is the bounds on the behavior of a solution to such 
a PDE. One classical result is that solutions to the heat equation have a Harnack- 
type inequality which says that the maximum temperature in a domain R can be 
known in terms of the minimum. That is 

max u(-, ti) < Cmin u(-, t 2 ). 
R R 

In this inequality, the constant C can be calculated explicitly based on the first 
time observed, t±, the second time observed, t 2 , and the dimension of the space in 
which the heat diffuses. A second property that is often of interest is how smooth 
or well-behaved the solutions are. One goal is to show that solutions are Holder 
continuous with respect to some metric. If x and y are any two points in the domain 
considered, then 

\u(x)-u(y)\ <C\x-y\ a 

where a is between zero and one. 

It is natural to extend the restrictions that one places on the equation to gain 
further understanding of the nature of the diffusion of heat. We may add the 
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complication of specifying the nature of the medium through which the heat diffuses. 
We do this with a matrix that models the homogeneity of the material. The new 
equation would be 



where a = is an n x n matrix of functions that model the material. This 

equation is well known when we assume the coefficients are smooth. This would be 
the easiest case. Even in the worst case, where the functions are merely measurable, 
much is known. 

A natural extension of the original problem would then be to take the derivatives 
not along the coordinates axes, but rather along vector fields. The vector fields 
could yield strange curves and even fractal-type behavior. In general, we consider 
less vector fields, m than the dimension of the space, n. But in order to do much 
with them we ask that they possess a certain property: that they can span the 
entire space being considered when taken in any combination. If X { is any of the 
m vector fields, then the equation takes the form 



Because the solutions cannot be assumed to be too smooth, we add another com- 
plication. We assume that the derivatives do not necessarily exist as functions but 
only in a weak sense. This further muddies the original problem by making it more 
difficult to take the two derivatives we have in the above equations. 

A final wrinkle on the original equation would be to consider not a matrix a^, 
but a non-linear function A in the equation. This gives us 



u t = div (aVw) 



rn 




m 



j 
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Non-linear behavior is often difficult to deal with in PDE. This last change would 
give us a much wider class of functions to consider. 

This paper proves the two properties considered above regarding the behavior 
of solutions. We show a Harnack inequality and Holder continuity of the solu- 
tions. We assume that the coefficients are only measurable, bounded functions and 
not smooth. We assume the derivatives are taken along vector fields that yield 
potentially strange geometry. 

This paper deals with two separate but related problems. First we will consider 
positive, weak solutions, u(x,t) G S 1,2 (R), to 



in a domain R = B cc (x, r) x (0, r) C R n+1 where the {X i }'^ =1 , m < n, are Horman- 
der vector fields. The matrix of coefficients {aij(x,t)} are bounded, measurable 
functions on R satisfying an ellipticity condition. We will call this the parabolic 
problem. For this problem we deal with two issues, establishing a Harnack-type 
inequality for positive, weak solutions and proving Holder continuity for those so- 
lutions. The problem is a classical one addressed by many authors. Jurgen Moser 
[25J dealt with weak solutions, u(x,t) G H /1,2 (Q), to 



in a unit cube, Q C M. n . In this paper, he introduced an iteration scheme that has 
formed the foundation of an approach to the issue of regularity of solutions adapted 
by many authors. His general approach was to define for some bounded domain 

D C R n , 






Using the facts that 



and 



maxn = M(oo, D) 



min u = M(— oo, D) 



D 

combined with Cacciopoli and Sobolev inequalities he was able to construct an 
iteration that accomplishes for some constant C > 



max M = M(oc,D) = lim (-L / u pk ) 1/pk < C (jL [ u ph ) 1/pk 
along with 

min M = M(-oo, J D)= lim [ u^Y^ > C f u^ k f hpk . 

d k^oo y \D\J D 1 y \D\J D ' 

These two inequalities hold for values of pk close to 0: — pk < < pk. Then these 
are linked with 

(JL f u p>*f pk < c (rL [ u - pk ) lhPh 

\D\ Jd \D\ J d 

for pk — > e. Moser later |26j addressed the continuity of solutions to a parabolic 

operator in W 1 using the same iteration. We will employ the Moser iteration here. 

The second problem we deal with is positive, weak solutions to 

m 
3=1 

in the same region R as above. Here, {Xj}™^, m < n, are smooth vector fields 
satisfying Hormander's condition for hypoellipticity. X* is the formal adjoint of Xj, 
and A(x, Vo«) is a measurable function satisfying certain requirements (see (j0.5p 
and following). This is a nonlinear generalization of the parabolic problem. It is 
useful to extend the Moser technique to many situations as authors have since his 
work in 1961. Trudinger [36] used the Moser method in 1967 to establish Harnack 
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type inequalities for quasilinear elliptic equations of the form 

div a(x, u, Vw) + b(x, u, Vu) = 

and 

div a(x, t, u, Vw) + b(x, t, u, Vw) — u t = 0. 

Then in 1995 Saloff-Coste, [32J used a sketch version of Moser's method to show 
that the parabolic Harnack inequality is equivalent to both the doubling property 
for measures and the Poincare inequality in the Riemannian manifold setting. This 
small list is far from exhaustive. 

We note that the following results are motivated by a direct proof, using the it- 
eration of Moser, of the Harnack inequality and Holder continuity of weak solutions 
to the two equations we examine. In the case of derivatives along Hormander vector 
fields, results are known in the case of certain elliptic operators, (see [7J) but not 
in the parabolic case. Also, results are not known for the case where the matrix of 
coefficients only measurable and L°° bounded. Results have been known for 

parabolic-type operators since 1958 with constant a^- coefficients (see |16j ). Then 
for measurable coefficients a parabolic and elliptic Harnack inequality was proven 
for non-divergence form operators (see [21]). This current work brings together 
the following: derivatives along Hormander vector fields, merely bounded and mea- 
surable coefficients, parabolic-type operators and Moser's method. To accomplish 
this new tools are needed. A Sobolev inequality for Hormander vector fields from 
Capogna, Danielli, Garofalo, (see [7]), a Poincare inequality for Hormander vectors 
fields from Jerison, (see [19]) and the existence of smooth test functions of compact 
support on the Carnot-Caratheodory balls with bounded horizontal gradient from 
Citti, Garofalo and Lanconelli, (see |10| ). The iteration scheme of Moser is nec- 
essarily different. We take into account the geometry of the Carnot-Caratheodory 
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balls that makes the integration tricky. These all result in explicit constants for 
both the Harnack inequality and Holder continuity in this geometry. 

0.2 Main results 

The main results of this paper are the following four theorems. For the parabolic 
case we have the 

Theorem 0.1. Harnack inequality for weak solutions 
Let u G S 1,2 (R) be a positive, weak solution to 

m 

(0.1) u t = -Y i X;(a ij (x,t)X j u) 

in R = {B cc (x,r) x (0, r)}. Also define R + = {B(x,r') x (r + ,r)} and R~ = 
{B(x,r') x (rf,r 2 -)} for < r> < r, < rf < r 2 < r+ < r. The {X,}™ 1; 
m < n, are smooth Hormander vector fields and the matrix of coefficients a,ij(x,t) 
are L°°(R) bounded and elliptic satisfying 

(0.2) A- 1 |e| 2 <a,,(x,t)^<A|e| 2 

for every (x,t) G R and for every £ G R™. Then there exists a constant C > 
depending on A > 0, n, and R such that 

(0.3) max-u <Cmin-u 

R- R+ 

where the maximum and minimum are the essential supremum and infimum over 
their respective domains. 

We use the notation for some function v, 

X*v = — > (tt~ v + — bA. (See section 0.5) 



1 dxi dxi 



The second is the 
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Theorem 0.2. Holder continuity of weak solutions with respect to the CC metric 

Let u(x,t) G S 1,2 (R) be a positive, weak solution to 

m 

u t = ~Y^ X*(a ij (x,t)X j u) 

in R as defined above. Then for every (x,t) and (y,s) G R = B(x,r) x (0, r), 
< r < 1 there exists a constant C depending on n, A, and R such that 

(0.4) \u(x,t)-u(y,s)\< Cd((x,t),(y,s)) a 

for some < a < 1 and d{[x, t), (0, 0)) = max {d cc (x, 0), v^RI}- 

Then for the nonlinear generalization we have the following 

Theorem 0.3. Harnack inequality for weak solutions to a nonlinear operator 
Let u G S 1,2 (R) be a positive, weak solution to 

m 

(0.5) u t = -Y,X;A j {x,V u) 

3=1 

in R= {B{x,r) x (0,r)}. Also define R + = {B(x,r') x (r+,r)} for < r' < r, 
< Tj - < tvT < t + < r. m < n, are smooth vector fields satisfying 

Hormander's condition for hypoellipticity . X* is the formal adjoint of Xy The 
measurable function A(x, Vqu) fulfills the following structural requirements: 



(S) 



Aj(x,V u) < C\V u\ 
\\V u\ 2 < Aj(x,V u)Xj < A|V m| 2 



for some constants C , A, A > 0. Then there exists a constant C depending on \, 
n, and R such that 



(0.6) 



max u < C min u 
R- R+ 
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where the max and min are the essential supremum and infimum over these do- 
mains. 



Last we have 

Theorem 0.4. Holder continuity of weak solutions to a nonlinear operator with 
respect to the CC metric 

Let u(x,t) G S 1,2 (R) be a positive, weak solution to 

m 

ut = -J2 X j A A x ^ou) 

3=1 

in R as defined above. Then for every (x,t) and (y,s) G R = B(x,r) x (0, r), 
< r < 1 there exists a constant C > depending on n, A, and R such that 

(0.7) \u(x,t)-u(y,s)\< Cd((x,t),(y,s)) a 

for some < a < 1 and t), (0, 0)) = max {d cc (x, 0), a/R[}- 

0.3 Examples of the nonlinear parabolic operator 

Let u G S 1,2 (R) be a positive, weak solution to 
(0.8) u t =-X;A j (x,V u) 

in R = {B{x,r) x (0,r)}, R C R n . Also define R + = {B{x,r') x (t+t)} for 
< r' < r, < rf < r 2 _ < r + < r. {Xj}™ 1; m < n, are smooth vector fields 
satisfying Hormander's condition for hypoellipticity. That is: 

Rank (Lie[X 1? X m ] )(x) —n 
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for any x G M. n . X* is the formal adjoint of Xj. The measurable function A(x, Vou) : 
R m — > IR m fulfills the following structural requirements: 




A,-(z,V u) <C\V u\ 
A|V n| 2 < Z7=iMx^ u)X jU < A\V u\ 2 



for some constants C, A, A > 0. 



i. First we point out that the parabolic version we prove in Theorem 10.11 is 
contained in this nonlinear example. Let Aj (x, Vo«) = a i j{x 1 t)X i u. Then 
we have 

u t = -) j X*A j (x, V u) = - X* (ajj (x, t)Xju) . 

Then we have the structural requirements (S) met in the inequality (j0.2p 
from the ellipticity and boundedness of dij(x,t): 

and noticing that 

aij(x,t)XiU < X\XiV,\. 

We can apply the results of Theorems 10.31 and 10.41 to achieve the same 
Harnack inequality and Holder continuity proved directly in the parabolic 
case. 

ii. Let 

Aj(x,V u)- 



Vl + lVowl 2 

with |Vo"u| < M in R for some constant M > 0. This is similar to the 
mean curvature for intrinsic minimal graphs studied by several authors: 
Cheng, Hwang, Malchiodi and Yang [9J, Manfredini [24j and Pauls [31J. 
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We consider weak solutions of the equation 



Ut = y Xj ( - ^ a ). 

The structural requirements are met in the following: 

A x,V w -V w = , = , ■V w< J -y- 
v/1 + |V u| 2 VI 

and 

A>,V n)-V n= - VoM -V 0M > |VoM|2 



1 + |V u| 2 ~ vTTm 2 



Applying the results from Theorem 10.31 we have a Harnack inequality for 
weak solutions to (10. 5p in R. We also are guaranteed Holder continuity of 
these solutions with respect to the CC metric, 
iii. Let 

with |Vqw| < M in R for some constant M > 0. That gives us 

U t ~ v V ,, , — 



Y / X*A j (x,X j u) = -J2x* T 



|V w| 2 

We can see that the structural conditions are satisfied in the following: 

a / „ x -XiW IVnwl .„ 

A(z,V w) = — , < - , ° 1 „ <|V w|, 



A(ar, Vo«) ■ Vo^ = Aj(x, Vqu) ■ XjU 



Vpw 
- |V w| 5 



V w < 



|V wp 



11 



and 



m 



V m 



|Vo^| 2 
1 + |M| 2 ' 




l + |V w| 2 



■ V m > 



Applying the results from Theorem 10.31 we have a Harnack inequality for 
weak solutions to (10. 5p in R. We also are guaranteed Holder continuity of 
these solutions with respect to the CC metric. 

0.4 Chronology of the elliptic and parabolic Harnack inequality 

This is a table showing some developments relevant to the current work in the 
study of elliptic and parabolic operators. We will look at different settings and with 
assumptions on the matrix of coefficients, A(x,t), having components ciij(x,t). 

1954 - J. Hadamard, [16J shows the first parabolic Harnack inequality for 
operators with constant coefficients in ~R n . Pini obtained the same result in 1954 
separately. 

1958 - J. Nash, [30J proves a Harnack inequality for an elliptic operator in M n 
using heat kernel estimates. This was for measurable and uniformly elliptic {(%}. 
Moser showed the same in 1964 using the iteration method that is employed in the 
current paper. 

1961 - J. Moser, [26j proves an elliptic Harnack inequality for solutions to 



and that the Harnack inequality implies Holder continuity. 

1964 - J. Moser, [25J shows a parabolic Harnack inequality using his iteration 
scheme and the John-Nirenberg theorem. This paper provides the basis for many 
authors ( see |12| . [36) . [32j ) as well as the current results. 
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1967 - Aronson and Serrin, [2] treat the second order quasilinear parabolic 
equation 

u t = div A(x, t, u, Vtt) = B(x, t, u, Vw) 

and were the first to show two-sided Gaussian bounds for divergence form uniformly 
elliptic operators in IR n . 

1967 - Trudinger, |36j deals with two kinds of equations 

div a(x, u, Vu) + b(x, u, Vw) = 

and 

div a(x, t, u, Vw) + b(x, t, u, Vm) — u t = 

and following the iteration used by Moser establishes a Harnack inequality for weak 
solutions to both. 

1968 - Trudinger, |37| shows point wise estimates for quasilinear parabolic 
equations. 

1969 - Bony, [3] considers degenerate elliptic partial differential equations of 
the type 

r 

X?u + Yu + cu = 

i=l 

where Xi, X T , Y are first order homogeneous differential operators. Various prop- 
erties of solutions are studied in connection with the Lie algebra generated by these 
operators. He proves a local form of the Harnack inequality. 

1971 - Moser, [27J shows a parabolic Harnack inequality using alternate meth- 
ods besides the John-Nirenberg theorem. The John-Nirenberg theorem is used in 
the current results in a form modified by Aimar pp. Moser uses an approach by 
Bombieri based on minimal surfaces. 
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1973 - Trudinger, |38j deals with operators of the type 

d df 
Lu = ——(aiAx)u Xi + aAx)u) + bi(x)u Xi + a{x)u = h / 

OXi v OXi 

where the measurable functions on a domain Vt C W 1 and / is an integrable 

function on Q. He proves the Harnack inequality for solutions and then shows 
Holder continuity. 

1980 - Krylov and Safonov, [21] prove a parabolic and elliptic Harnack 
inequality for non- divergence form operators 



Lu = — ciijdid, 

i,3=± 



with dij measurable, symmetric and uniformly elliptic. 

1983 - Pranchi and Lanconelli, |12| were the first to apply a Moser technique 
to obtain Holder continuity of solutions for non-uniformly elliptic operators with 
measurable coefficients. 

1986 - Jerison and Sanchez-Calle, |20| treat degenerate elliptic operators 
with smooth coefficients. We give their main result. 

Theorem 0.5. Let Q be an open subset opR n . Let 

l=x:(ia(-))^(^k/) 

i,j=l % x i 

where the coefficients and k are C°°(Q), k is positive and the matrix is 
symmetric positive semi-definite for every x £ Q. Also suppose that L satisfies a 
sub-elliptic estimate: There is a constant C and a number e such that allu G C^°(Q) 
satisfy 

\\u\\ 2 e < C{\\Lu\\ + \\ U \\). 

Here 

|2u jAl/2 



( / i«(or(i + iei 2 )*de) : 
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and 

dfi(x) = h(x)dx. 

Then for ( — L)u(t, x) = we have for the heat kernel h(t, x, y) the following: 

h(t } x } y) < A^{B{x,t l/2 )Y 1 e~ d{x ' y)2ht 

h(t,x,y) > AV(5(x,t 1/2 ))~V d( ^ )2/ ^ 
\d*X h ---X jp h(t,x,y)\ < C^r^^Bix^/ 2 ))-^-^ 2 ^ 
For some positive constants A, A', 7, 7', C Pi k and every x,y £ Q < t < 1. 

1988 - Kusuoka and Stroock, [22J use probability methods for heat kernel 
estimates on Hormander vector fields. We give their main result. 

Theorem 0.6. Let aij G C 2 be a symmetric, non-negative definite matrix-valued 
function. Define 

N 

Lu{x) = ^2[d Xi (a i:j d x .u))(x). 

Let P(t, x, •) be the unique transition probability function on W N such that the as- 
sociated Markov semigroup {P t : t > 0} satisfies 

P t( j>{x) - <p(x) = [ [P s L<P](x)ds 
Jo 

for all cf) G Cq°(M. n ). Let {E\ : A G [0,oo)} denote the resolution of the identity 
determined by {P t : t > 0} and set A = Jj 0oo N \dE\. Let the Dirichlet form be given 
by 

S(fJ)= [ \d(E) 
J [0,00) 
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Assume there exist A G (0, oo), v G (0, oo) and 5 G (0, oo) such that 



<A(S(f,f) + 5\\f\\l)\\f\\r 

for f G L 2 (M. N ); or equivalently that there is a B G (0, oo) such that 

\\Pt\\i^oo<Be H /t»/ 2 

for t > 0. Then P(t, x, dy) = p(t, x, y)dy and there is a C G (0, oo) depending only 
on v such that for each p G (0, 1] and all (t, x) G (0, oo) x M. N : 

pit, x, •) < C(A/pty/ 2 e pst exp[-D(x, -) 2 /(l + p)t] a.e. 

1992 - Varopolous, [41J shows various parabolic Harnack inequalities using 
heat kernel estimates. He does not use the approach taken by Moser. 

Theorem 0.7. Let {Xi}^ =1 be a Hormander system of vector fields on the manifold 
V , Y a C°° vector field on V , and a, a C°° function on V. Let D be the differential 
operator 

i=l 

defined on K x V. Let I = (a, (3) C M., Q and open relatively compact subset of V , 
K a compact subset of Q, t\, t 2 , such that a < t\ < t 2 < (3, J G I{N), m G N. 
Then there exists a constant C such that every positive solution u of Du = in 
I x Q satisfies 

SU P l(|:) m (|-) J ^1^)1 ^ Cm f u ^ x ) 

x eK ot OX xG k 
with spatial derivatives taken with respect to some fixed system of local coordinates. 



1992 - Grigory'an, [15J uses heat kernel estimates to prove the parabolic 
Harnack inequality and does not rely on Moser's method. He shows the Harnack 
inequality is implied by the doubling property and the Poincare inequality. He 
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also shows a Harnack inequality without using the doubling property and Poincare 
inequality. The behavior of the Green function G(x, y, t) of the Cauchy problem for 
the heat equation on a connected, noncompact, complete Riemannian manifold is 
investigated. For manifolds with boundary it is assumed that the Green function 
satisfies a Neumann condition on the boundary. His results are similar to Saloff- 
Coste [32j but arrives at them with different methods. 



Theorem 0.8. Let M be a geodesically complete, noncompact, smooth, connected 
Riemannian manifold of dimension n. Let A be the Laplace operator on M as- 
sociated with the Riemannian metric g. Suppose M satisfies the following two 
hypotheses: 

1. For any concentric geodesic balls B(x,R), B(x,2R) of radii R and 2R 



Where the constant A is the same for all balls. 

2. For some constant N > 1 and a > 0, for any geodesic ball B(x,R) C M 
and for any smooth function f defined in the ball B(x, NR) the following inequality 
holds: 



Then the uniform parabolic Harnack inequality is valid with the constant P depend- 
ing on A, a, N . 

Fix a point x G M . Let u(x,t) be a positive solution of (d t — A)u = in Q 8 r = 
B(x,R) x (0,8-R 2 ) which is smooth in Qsr and satisfies the Neumann condition 
for x G dM (if dM is non-empty). Set Q = Br x (3R 2 ,4R 2 ), and suppose that 
sup u = 1 . Then u(x, 64R 2 ) > 7, 7 = j(A, a, N) . 



HB(x,2R) < Afj,(x,R) 




Q 
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1993 - Citti, Garofalo, and Lanconelli, [10J establish a uniform Harnack 
inequality for weak solutions u of 



where V is a measurable function belonging to an analogue to the local Kato class 
for the Xj. They also introduce a very important cutoff function for CC balls which 
we make use of in the current paper. We state this lemma for completeness. 

Lemma 0.9. Existence of cut-off functions 

For every ball B(x, t) C R, and every s,t > with s < t, there exists a function 
ip G C$°(B(x,t)) such that < *0 < 1 in B(x,t), ip = 1 in B(x,s) and \V ^\ < 
C/(t — s), where C is a positive constant which is independent of t and s. 

1993 - Capogna, Danielli and Garofalo, [7] prove a version of the Sobolev 
embedding theorem for functions in Sq ,p using the fractional integral operators used 
by Gilbarg and Trudinger [14] . With this result they prove a Harnack inequality 
for solutions to a subelliptic operator and Holder continuity. Their main result: 

Theorem 0.10. Let U C R n be a bounded open set with relative homogeneous 
dimensionQ > 0. Let {X i }'^ =1 be C°° vector fields satisfying Hormander's condition 
for hypoellipticity . Assume 1 < p < Q. Let u G S]f c be a non-negative solution to 
the equation 



Then there exist C > and R > such that for any Br = B(x, R) with B(x, 4R) C 
U, and R < Rq, 



rn 



u = 



in 




ess sup u < C(ess inf u + K(R)) 

Br B r 
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1993 - Varopolous, [41J shows that the upper and lower Gaussian bounds on 
the fundamental solution are implied by the parabolic Harnack principle. This is 
done in the setting of Hormander vector fields. 

1995 - Saloff-Coste, |32j shows the parabolic Harnack inequality is equivalent 
to the doubling property and the Poincare inequality. This is done in Riemann- 
ian manifold setting. He assumes the a^- coefficients are bounded, measurable, 
symmetric and elliptic. His main result: 

Theorem 0.11. Let M be a C°° connected manifold. Let L be a second order 
differential operator with real C°° coefficients on M and no zero order term. Use 
as distance p(x,y) = sup{\f(x) — f(y)\,f G C°°(M), |V/| < 1}. Consider the two 
properties, the doubling property and Poincare inequality: 

(0.9) K B ( X > 2r )) < CiK B ( x i r )) 

(0.10) f \f-f B \ 2 dp<C 2 r 2 [ |V/| 2 rf/i 

J B(x,r) JB(x,2r) 

for0<r<r ,xe M, f E C°°(M) and f B = js\ f B f 
Then the following two statements are equivalent 

A. The properties W.9\) and A0.10\) above hold for some r . 

B. There exists r\ > and a constant C depending only on the parameters 
0<e<r]<5<l such that for any x G M, any real s, and any < r < T\, any 
non-negative solution u of (dt + L)u = in Q = (s — r 2 , s) x B(x, r) satisfies 

(0.11) sup u < Cinf u 

Q- Q+ 

where 

Q_ = [s — 5r 2 , s — r]r 2 } x B(x, Sr) 
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and 

Q + = [s — er 2 , s) x B(x, Sr) 

It should be noted the similarities and differences between this paper's current 
results and those of Saloff-Coste. The current results are contained in those of 
Saloff-Coste. We examine a parabolic-type operators that satisfy properties (10. 9p 
and ( IQ.lOp above. So by his results, the Harnack inequality holds in some region. 
He does not employ a direct proof which we do here. Also, using the geometry 
created by the Hormander vector fields, we calculate the constants for the Harnack 
inequality and Holder continuity explicitly. A further advantage of the direct proof 
is to allow for extension into the generalized nonlinear case as we do here as well 
as other cases. 

2002 - Uguzzoni, Bonfiglioli, and Lanconelli, [39J use heat kernel esti- 
mates with Gaussian bounds to prove the parabolic Harnack inequality on Carnot 
groups. 

2004 - Uguzzoni, Bonfiglioli, and Lanconelli, |40j construct the funda- 
mental solutions for non- divergence form operators where the Xi are Hormander 
vector fields generating a stratified group G and positive-definite matrix 

with Holder continuous entries. They also provide Gaussian estimates of T and its 
derivatives and some results for the relevant Cauchy problem. 

2006 - M. Bramanti et al., [4J prove the existence of a fundamental solution 
for a class of Hormander heat-type operators. For this solution and its derivatives 
they obtain sharp Gaussian bounds that allow to prove an invariant Harnack in- 
equality. The matrix is real symmetric, uniformly positive definite and having 
Holder continuous entries. 

A list of some developments relevant to the current work of those in the Hor- 
mander vector field setting. Chronologically: 

1969 - Bony showed local Harnack inequalities for degenerate elliptic operators. 
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1988 - Kusuoka and Stroock use probability methods for heat kernel estimates. 

1992 - Grigory'an shows a Harnack inequality without using the doubling prop- 
erty of measures and the Poincare inequality. 

1993 - Capogna, Danielli and Garofalo prove a Harnack inequality and Holder 
continuity for solutions to a subelliptic operator. 

1993 - Varopolous shows that Gaussian bounds on the heat kernel are implied 
by the parabolic Harnack principle. 

1995 - Saloff-Coste shows the parabolic Harnack inequality is equivalent to the 
doubling property and the Poincare inequality. 

2002 - Uguzzoni uses heat kernel estimates with Gaussian bounds to prove a 
parabolic Harnack inequality. 

2004 - Uguzzoni constructs the heat kernel for the non-divergence form opera- 
tors. 

2006 - Bramanti uses sharp Gaussian bounds to prove a Harnack inequality. 
0.5 Hormander Vector Fields 

We will use vector fields of the form 

OXi 

where the b{ G C°°(lR n ). We will consider solutions to 

m 

u t = - ^ X*(a ij (x,t)X j u) 

and 

m 

along smooth vector fields {Xj}™ 1; m < n, satisfying Hormander's |18| condition: 

rank[ Lie[Xx, X m ]j =n 
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at every x G M n . 

We may take the derivative along a vector field. If we consider a simple example: 



then we have for some function /: 

OX i 

We also will make use of the formal adjoint of any JQ, X*. The formal adjoint of 
a vector field Xi is the operator X* such that for functions /, g G C^°(R) 

[ X i f-gdxdt= [ f-X*gdxdt. 
Jr Jr 

We compute: 

/ X i f-gdxdt= / yV& j75 -/).<7<fa;df= / ^(W) • (6^) ttedt 

</R </R J=1 i=1 ° X 3 

And integrating by parts gives us: 

= - 1 J • ^ = - / r / . + ^) *** 

= [ f-X*gdxdt 
Jr 

So that X*g takes the form (-l)(^g + J^fy) 
In this paper, we use equations of the form 

m 

u t = -J2 X j A A x iV u) 

3=1 

So for G C£°(R) we can say 

m 
3=1 
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And integrating over R 

„ m „ 

/ 4>u t dxdt = — / Xj(f)Aj(x,Vou) dxdt. 
Jr j=1 Jr 

We use throughout this paper the horizontal gradient of a function /. 

Definition 0.12. Horizontal gradient 

The horizontal gradient of a function f is defined to be: 

V f = X l fX 1 + X 2 fX 2 + ■■■ + X m fX m . 

And the size of the horizontal gradient is 

m 

iv /i = (B^) 2 ) 1/2 - 

1=1 

We note that for the the horizontal gradient will not in general equal the 

Euclidean spatial gradient in R™, V/. The vector fields {X i }™ 1 give us a natural 
distance. 

0.6 Carnot Caratheodory Distance and Balls 

The vector fields {Xi}^ give us a control distance. Consider a piecewise C 1 
curve 7 : [0,T] ->• lR n . Let x, y E R n . For 5 > we define the class C(5) of 
absolutely continuous paths 7 : [0, 1] — > 1R 3 with endpoints 7(0) = x and 7(1) = y, 
so that 

m 

1=1 

and 

m 

^a,(t) 2 <5 2 

i=l 

for a.e. t G [0, 1]. Such paths are called horizontal. 



Definition 0.13. Carnot- Caratheodory metric 



We define the Carnot-Caratheodory (CC) metric to be 
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(0.12) 



dcc(x, y) 



inf{5 such that C(6) = 0}. 



A dual formulation is 



d cc (x, y) = inf {T : 3 7 : [0, T] R n , 7 (0) = x, 7 (T) = y,J2 ^ 1 a ' e ->- 



We will refer to the CC distance as d and d cc . With the CC distance we can define 
balls in R n . We set 



the ball centered at x of radius r. We use as the volume of the ball \B(x,r)\ = 
Ib(x r) d- x - ^ is worth noting that the geometry created by the vector fields, {Xj}^ 
is different from the Euclidean space, MJ 1 . One major issue we face in obtaining 
estimates on the behavior of functions is the geometry of the Carnot-Caratheodory 
balls. In the Euclidean setting the ratio of the volume of two balls is proportional 
to their radii. In the Carnot-Caratheodory setting the situation is different. The 
volume of two balls is related in the following way from a proposition of Nagel, 
Stein and Wainger [29J. 

Proposition 0.14. Given a bounded set U C W l there exist Q > n, Rq > 0, 
Cb > 0, such that for every x G U , R < Rq and < t < 1 



which we will make use of in our results. Test functions play an integral role in 
these results. Since the Sobolev theorem we use is assumed to be for functions of 
compact support on a Carnot-Caratheordory ball and our solutions are not sup- 
ported in these domains, we work around the problem with certain test functions. 



B r = B(x,r) = B cc (x,r) = {y : d cc (x,y) < r} 



(0.13) 



B(x,tR)\ > C B t Q \B(x,R)\ 
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Moser used such functions in the Euclidean case [25J, |26j and they are standard 
in applications of this type. But the geometry of the Carnot-Caratheordory balls 
is much different. In the Euclidean context, balls in n-dimensions are very well 
behaved, having symmetry as well as level sets which are also balls. This is not 
the case in the present context. The Carnot-Caratheordory balls can be badly 
behaved, having cusps and the level sets could potentially behave just as badly. 
When we use the test functions, we set two radii and construct the function to be 
= 1 on the smaller of the concentric balls. Then the function decreases to with 
bounded gradient until it is <ft = outside of the larger radius. This fact has been 
an obstacle until 1993 when Citti, Garofalo, and Lanconelli jlOJ constructed such 
test functions. Also Nagel and Stein [28J in 2001 revisited the test functions along 
with showing the existence of smooth metric equivalent to the control metric over 
Hormander vector fields. 

There is another difficulty present in the use of these test functions in the 
Carnot-Caratheodory setting. In the construction, the two radii may not be ar- 
bitrarily close to each other. We may not merely choose two radii r > r' > 
and construct a useful test function. The radii must be bounded away from each 
other. We recall Nagel and Stein [Lemma 3.1.1] [28J that there exist two constants 
< Ci < 1 < C2 such that = 1 on x < C\T and = on x > c^t. Here c\ and c<i 
come about in the construction of level sets of the CC balls. 

The facts of this construction are worth mentioning even though the details of 
these restrictions are not obvious when we use the test functions. In constructing 
the test functions Citti, Garofalo and Lanconelli construct a quasi- distance D(x, y) 
that satisfies 

cid cc {x,y) < D(x,y) < c 2 d cc (x,y) 
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for every x,y £ B(x,r). Now construct a new ball B(x,R) = {y | D(x,y) < R}. 
This gives us 

Br C B(x, R) C Br 

which says that we can create a B ball whose boundary is contained in a CC ball 
of a larger radius and containing a CC ball of a smaller radius. In order for us to 
use this fact we must construct nested balls. In the proof we use radii r and r' so 
let those constants be given. This gives us two inequalities based on 

dd cc (x,y) < D(x,y) < c 2 d cc (x,y) 

which are 

C]T < r < c 2 r 

And the second is 

C\r < r < c 2 r . 

So we need balls that satisfy 

B(x,r') C B(x,r'c 2 ) C B(x,rci) C B(x,r). 

This requires our constants to satisfy r'c 2 < rc\. Recalling that < C\ < 1 < c 2 
forces r and r' to be bounded away from each other and not arbitrarily close. 

0.7 Weak solutions 

We will consider positive, weak solutions to 

m 

u t = - ^2 X*(a ij (x,t)X j u) 



and 



3=1 
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in a region R C R Tt+1 . In order to create the context for these solutions we will 
begin with definitions. 

Definition 0.15. L 2 (R) 

L 2 (R) is defined to be the space of all functions such that 



i/iu*(*) = ( / m 2 cM) 1/2 <oo 

Jr 



Definition 0.16. L 2 oc (R) 

Lf oc (R) is defined to be the space of all functions such that 

II/IIl»(*> = (/ Iff dxdt) 1 / 2 <oo 

J K 

for all compact subsets K of R. 

Throughout the paper we use the notion of weak derivatives. We define 

Definition 0.17. Weak derivative 

If u, v G L 2 (R), then v is the weak derivative of u if 



I uX*4> dxdt = — V(f> dxdt 
Jr Jr 



for all (f) e C%°(R) 



Definition 0.18. S^R) 

We define S 1,2 (R) to be the set of functions u : R — > R such that u, u t are in L 2 (R) 
and all of the horizontal derivatives of u are in L 2 (R). 

The S 1,2 (R) norm is then given by 

rn 

||ti||sl,2 = ||u||l 2 + \\ u t\\ + ||^w||l2. 

i=l 

We note that S 1,2 (R) is a Banach space under this norm. 
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In this paper we consider weak solutions. 

Definition 0.19. Weak solution to $0.1\) 
A weak solution to 

m 

u t = - 22 X*(a ij (x,t)X j u) 

is defined as a function for which the first derivatives Ut,X\u, ...,X m u are square 
integrable in R and satisfy 

„ „ m 

/ / (j>u t + V"^ (aij(x,t)XjUj X; L (f) dxdt = 
JJr i,j=i 

for every 4>(x,t) G C°° which has compact support in the x variable for every fixed 
t. 

We use an analogous definition for weak solutions to ( 10. 5 j) 

Definition 0.20. Weak solution to $0.5\) 
A weak solution to 

m 

i=i 

is defined as a function for which the first derivatives Ut,XiU, ...,X m u are square 
integrable in R and satisfy 

/ / (j)u t + (Aj(x, Vow)) Xj(f) dxdt = 

JJr j= i 

for every (j)(x, t) G C°° which has compact support in the x variable for every fixed 
t. 



CHAPTER 1 
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Parabolic Harnack Inequality 



1.1 Flow of the proof 

The flow of the proof will be as follows. We will first prove a Cacciopoli inequality 
for weak solutions v to ( 10. ip . We will show that these hold for various powers of the 
function. This will give us two inequalities for bounds on both the square of the 
horizontal gradient of v and on the maximum of v 2 over all t from above in terms 
of v 2 . We will the use a Sobolev embedding theorem by Capogna, Danielli and 
Garofalo [6j to allow for the estimate of higher powers of v in terms of v 2 . We will 
combine this with the bounds from the Cacciopoli inequality to form an inequality 
that we may iterate to achieve both the minimum and maximum of our function 
over specific domains. The iteration will work for integrals of v p and v~ p for p close 
to 0. We will bridge the gap between these two values with the results of Aimar [TJ 
and his John-Nirenberg style theorem for spaces of homogeneous type. In order to 
show that his results work, we need to show that log(v) satisfies a certain bounded 
mean oscillation (BMO) condition. With this result, we will finish the proof of the 
parabolic Harnack inequality. Then, patterning our proof after Moser |25j we will 
show that this Harnack inequality implies Holder continuity of solutions in a unit 
rectangle with respect to the CC metric. 



In order to begin this proof we will construct a Cacciopoli inequality relating 
the value of J RI v 2 over R' C R C W n+1 to J R \V v\ 2 and J R v 2 . 

Theorem 1.1. Let v e S ,1,2 (i?) be a positive, weak subsolution to 



1.2 Cacciopoli Inequality 



711 



(1.1) 
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in R = {B(x,r) x (0,r)}. Also define R! = {B(x,r') x (0, r')} for < r' < r, 
< t' < t. Let {Xj}™ 1; m < n, be smooth Hormander vector fields satisfying 
Hormander' s condition for hypoellipticity. Let the matrix of coefficients aij(x,t) be 
L°° bounded and uniformly elliptic. That is 



for every (x,t) G R and for every £ G M ra . Then there exist constants C\, Ci 
depending on \, n, and R such that 



Where the max is the essential supremum over this interval. 

Consider the function (f)(x,t) = v(x,t)tp 2 (x,t). We will make the support of 
ijj(x,t) compact in the x variable for any fixed t and differentiable. Notice this 
also makes > 0. In order to accomplish the construction of we will recall the 
existence of test functions in Citti, Garofalo, and Lanconelli [10J 

Lemma 1.2. Existence of Carnot-Caratheodory test functions. 

There exists ro > such that given a metric ball B(x, r) CC R, with r < r^ and 
< r' < r, there exists a function ip G C^°(B(x,r)) such that 0<ip<l,ip = l in 
B(x,r') and IVq^I < Here, Cue > is a constant independent of r and r' . 



(1.2) 




(1.3) 




Define ip(x,t) := ipi(t)ip2{%) m the following way: 



(1.4) Mt) 



1 if r' < t < r 
if t < 



1 if d cc (x, 0) < r' 
if <i cc (^, 0) > r 

And ipi(t) is linearly interpolated from to r'. Since v is a subsolution to ( 10. ip 
i? we have 

(1.5) J J (pv t dxdt + yy Xi(f) (aij(x,t)Xjv) dxdt < 



(1.6) y y i) 2 vv t dxdt + J J x i( v ip 2 ) a ij (x,t)X j v dxdt < 0. 

We calculate Xj(t>-?/> 2 ) = XiVip 2 + 2if)Xiipv and substituting 



(1.7) 



/ / ip 2 vv t dxdt + / / XiVip 2 dij(x,t)XjV dxdt 
JJr u=x JJ R 



t)Xjvdxdt < 0. 

Notice that i?/> 2 (i> 2 ) t = tp 2 vvt which gives 



(I- 



J J -ip 2 (v 2 ) t dxdt + J J XiVip 2 aij(x,t)XjV dxdt 



< —2 ^ ^ vipXiifj aij(x,t)XjV dxdt. 
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Now consider the RHS of the equation 



'1.9) vip Xjip a i j(x,t) X 



j v 



and the Schwarz inequality: 

m 

(1.10) \ y ^wl)X t ^a ij {x,t)X^\ < (i> 2 (V v,aV Q v)(Voi>,aV ip)v 2 ) 1 *. 
And using Young's inequality we get: 

(1.11) < ^(V v,aV v)iJ 2 + v 2 (Voij,aVo^). 
Now substitute back into the RHS of ( 11.81) to get: 

l \v 2 ) 

R & J J R 



1.12) / / -ip (v 2 ) t dxdt + 1 1 (V v 1 aV v)ip 2 dxdt 



(1.13) ~ 2 // Y^ V ° W ' aVo ^ rf:Erft + t;2 ( V oV' ) aVoV')- 

Now gathering the ifj 2 (VoV, clVqV ) terms onto the LHS of (jl.8p we get: 



1.14) 1 1 ^ 2 (v 2 ) t dxdt+ J J ^{V v,aV v)ip 2 dxdt 



< J J v 2 (V ^,aV ^)dxdt. 
Use the ellipticity and boundedness of (x, t) to get 

(1.15) J J Ujj 2 (v 2 ) t dxdt + JJ ^\V v\ 2 ip 2 dxdt < 2\JJ v 2 \V ip\ 2 dxdt 

Notice that ^(i/j 2 v 2 ) t = 'ip'iptv 2 + vv t ip 2 . And use this fact by adding v 2 i/jijjt to both 
sides: 
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(1.16) J J ]^ 2 {v 2 ) t dxdt + v 2 ^ t dxdt + J J -^\V v\ 2 ip 2 dxdt 

< 2\JJ v 2 (\V ^\ 2 + \^t\)dxdt 

This gives 

(1.17) J J ^ 2 v 2 ) t dxdt + JJ l.\V v\ 2 ^ 2 dxdt 



< 



2\JJ v 2 (\V ^\ 2 + \Wt\)dxdt. 



Now integrate this over R a = {B(x,r) xO<t<a<r} and drop the second 
integrand on the LHS. We note that 



JJ R ^l v o^lV«>o. 



We choose a so that 



(1.18) / v 2 \t= a dx > — max / v 2 dx. 

JB(x,r') 2 t€(0,<7) J B (x,r') 

This gives us: 

(1.19) If hv 2 tp 2 ) t dxdt<2\ II v 2 (\V Q ip\ 2 + dxdt. 
Now consider the 

(1.20) max / v 2 dx < 2/ v 2 \ t=a dx 

te(0,a) J B ( x y) JB(x,r>) 

which we constructed above, — 1 on R' and we fix t, so 



max / v 2 dx < 2 / v 2 \ t=a dx 

te(0,<r) J B (x,r') JB(x,r>) 
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< 2 / v 2 ^\ =u dx 

>B(x,r') 



< 2 / v 2 ^\ =u dx 

JB(x,r) 



And since ip vanishes at t = 



= 2 / v z i) z \t= a - v 2 ip 2 \ t=0 dx = 2 // {v 2 i) z ) t dxdt 

JB(x,r) JJRcr 

Now recall (I1.19P integrated over R a . We will combine that with the ^ factor as 
well as the fact that the above is four times the (j!.19p inequality. Putting these 
together gives: 

(1.21) max / v 2 dx < 8A // v 2 (\ V ip\ 2 + \ipip t \) dxdt. 

te(0,ff) jB{x,r') JjR a 

Then use the non-negativity of the RHS and take the integral over the larger rec- 
tangle R for 

(1.22) max / v 2 dx < SX v 2 (\V^\ 2 +\4>i) t \) dxdt. 

*6(0,<r) J B(x,r') JJr 

In creating the test function according to Lemma 11.21 above, we use the result 
on the bounds of the horizontal gradient of test functions in Citti, Garofalo, and 
Lanconelli jlOj in order to attain a constant independent of ip. Specifically, from 
Lemma [1.21 we use IVo^I < fr^"- This allows us to say 

(1.23) |Vo^| 2 + \m < ( J ^ f + -^-) 
Which gives us 

(1.24) max / v 2 dx < 8\( , Chg xo + — - — ) [[ v 2 dxdt. 

te(0,e)J RI V( r _ r /)2 t t J J ji 
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This proves (II .3p . Now look at the second integrand on LHS of (j!.17p . 
(1.25) jJ^^\V v\ 2 ip 2 dxdt < JJ^^\V v\ 2 dxdt 



R' 



(1.26) <2A // v 2 (\V ip\ * + \iH> t \)dxdt 

Which gives us 



1.27) ll \Vov\^ 2 dxdt<AX 2 JJ v 2 (\V ip\ 2 + \ipip t \)dxdt 



1.28) < 8A 2 ( , ^ 9 + - J—) / / v 2 dxdt 

v ' " \r-r') 2 t-t' J JJ r 



And the integral on the RHS is taken over R since it is non-negative. This gives 



1.29) UjV v\ 2 dxdt <8X 2 (-^y 2+: ^) /I v 2 dxdt 



with -0 = 1 on R' . 

1.3 Cacciopoli for other values of p 

Above we proved a Cacciopoli type inequality, (Theorem II. ip for positive, weak 
solutions to 

rn 

u t = - ^ X*(a ij (x,t)X j u) 

in R C IR n+1 for smooth Hormander vector fields. Now we need to show that this 
result holds for powers of u. We will consider a function v = u p for p > |. We also 
need similar inequalities for p < and < p < 1. This, along with the Sobolev 
style embedding inequality will provide the basis of the Moser-type iteration. 



Proposition 1.3. (Cacciopoli for powers of u) 
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Under the assumptions of Theorem \l.l\ above, the results are true for powers of 
v := u p , for p > \, 



1.30) J J \V (u p )\ 2 dxdt <d J J {vPfdxdt 



(1.31) max / {vPfdxdt < C 2 / / (u p ) 2 dxdt 

*6[0,<t] JB(x,r') JjR 

If u is a weak supersolution to W.l\) and p < then we have a similar result with 
the inequality reversed: 

(1.32) J J \V (u p )\ 2 dxdt >d j j '{vPfdxdt 

(1.33) max / (u p ) 2 dxdt > C 2 1 1 (u p ) 2 dxdt 

te[0,o] J B (x,r>) JJR 

Proof: 

Let v := u p be a positive subsolution to 

rn 

ut = - ^2 x *( a ij( x ^) x j u ) 

in the rectangle R as above and = ip 2 vP. Then we have 

(1.34) f ! <pv t dxdt+ ^2 ff X i (p(a ij (x,t)X j v)dxdt < 



R i,j=l JjR 



Calculate 



v t = (u p ) t = pu p 1 u t 



and 

Xi(<j>) = Xi(^ 2 u p ) = pu p ~ x i) 2 X iU + 2X i ipv p ip. 
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We substitute the definitions of and v. 

J j i?u v V u v - x u t dxdt 

+ V // {pvP-^XiU + 2X i ilJu p ^)a ij (x,t)(pu p ~ 1 X j u) dxdt < 



i) 2 u 2p ~ x u t + p 2 u 2p - 2 ^ 2 (V u, aV u)dxdt 

R 



+ jj 2pil)u 2p - l (V ^,aV u}dxdt < 0. 
We recognize that 2pu 2p ~ 1 u t = \ip 2 (u 2p ) t : 

J j ]^p 2 {u 2p ) t dxdt + p 2 u 2p - 2 ij 2 {V u, aV u) 

+ JJ 2pipu 2p ~ 1 (V 4>,aV u) dxdt < 
Examining the RHS we have by the Schwarz and Young's inequalities: 

(1.35) -2p^u 2p - 1 (V ^,aV u} < ftpim 21 " 1 (V ip , aV u}\ 

(1.36) < 2p((V ^, aV^)u 2p ^ 2 {V u, aV u}u 2p ~ 2 )^ 

(1.37) < 2p((V ^, aV ij)u 2p + -(V «, aV u)^fj 2 u 2p - 2 ). 
Substituting back in and combining like terms we have 

J j 4> 2 {u 2p ) t dxdt + JJ (p 2 - ^)u 2p - 2 4> 2 (V u, aV u) dxdt 
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< J J 2p(V ^, a\7 ip)u 2p dxdt 

At this point, observe v = u p , Vof = pvF^VoU and use both boundedness and 
ellipticity of dij(x,t) 

'.2/„.2\ J„J+ , / / 1 2 P ~ 1/.2IV7 „.|2 J„J + / / / r,„Mv7 „/.|2„,2 



(1.38) yy ^ 2 {v 2 ) t dxdt + J J - 1 ip 2 \V v\ 2 dxdt < J J 2p\\V ip\ 2 v 2 dxdt 

At which point, we have the same inequality as ( 11. 17ft and the proof follows as 
above. For the case p < po < we let v = uPl 2 be a supersolution to (10. ip in R. 
We then have for = v ip 2 the following: 

(1.39) v t = (u p/2 ) = V -u p ^u u 

(1.40) X t v = X t {u p ' 2 ) = V -u p l 2 ' x X lU 
and 



1.41) X i( f) = Xiivtp 2 ) = Xi(u p/2 ij 2 ) = ^u p/2 - l XM 2 + 2ij;Xrfm* 



) 2 ) = X.(n,P/ 2 nl, 2 \ = P-n.P/^Y-nnl, 2 -4- OnhY .<,hn,P/ 2 

And we substitute these into 



(1.42) J J (j)v t dxdt+^r J J X i 4 l (a ij (x,t)X j v)dxdt>0. 



to get 



+ 11 [^vr^XiUifi 2 + 2i)X i i)u p/2 } aij (x, t^vP'^Xiu) dxdt > 0. 
J jr 2 2 
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Since 



1.43) i(A = w*=(« p/a )(|« ?/a " 1 «*) 



we can arrive at 



-(v 2 ) t i) 2 + iJ 2 (V v, aV v) dxdt 
R 2 



> -2 wf>(V iJ),aVov) dxdt 



R 

At this point, we have inequality ( ll.8[) and we may proceed in the proof exactly as 
above. The case for < p < 1 is modified in the following way: we must reverse 
t and use v := u p (x, —t). Then our test function must be = u p ip 2 . The proof 
proceeds in the same way except there is a factor of j^r- in front of the IVo^] 2 
term. That is: 



1.44) / / ^ 2 {v 2 ) t dxdt+ 1 1 j 2p ^ 2 1 i) 2 \V qv\ 2 dxdt < / / 2p\\ V ip\ 2 v 2 dxdt 



1.4 Sobolev Theorem and Iteration 

Above we proved the following estimates, (Theorem II. ip : 

J J \V v\ 2 dxdt < C iJJ v 2 dxdt 

max / v 2 dx < C2 v 2 dxdt 

for positive, weak solutions to (10. ip in R. Then we showed that these held for 
powers of u, v := u v . Now we will combine these with the Sobolev embedding 
theorem to begin the Moser iteration. We introduce the 

Theorem 1.4. (Sobolev Embedding Theorem) Let U <zM. n be a bounded open set, 
then there exists a Cs > 0, r > such that for any x G U , B r = B cc (x,r), with 
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r < r o> here 1 < k < 2 



(1.45) 



I -^r | J B r \^r\ J B r 



for any u G Sr\' 2 (B r ) 



Squaring both sides, we will use this formula in the form 
(1.46) / \u\ 2k dx^ < C s (-L [ \V u\ 2 dx) 

I -°r | J B r I -"r I JB r 

In order to use the Sobolev embedding theorem in the setting where our function 
v is not compactly supported but merely defined on the ball B(x,r) we will recall 
the cut off function if)(x,t) defined above: 



(1.47) Mt) 



1 if t' < t < r 
if t < 



1 if d cc (x, 0) < r' 
if d cc (x, 0) > r 

and as before we have a bound on the horizontal gradient. So we set w = v ■ ip G 
Sr) ,2 (B(x,r)) and notice that on B(x,r') we have: 

w(x, t) = v(x, t) • ip(x, t) = v(x, t) ■ ip\{t) ■ ip2(x) = v(x, t) ■ 1 = v(x, t). 
Also we have 



V w = Vq{v ■ ip) = V vip + Voipv = Vqv>P + • v — Vovip = V f • 
Then on B(x,r) we have: 



40 



w(x, t) = v(x, t) ■ if)(x, t). 

Also we have 

V w = Vo(f • i>) = Vovip + Voi>v 
Now we calculate for the function w above: 



J- / \v\ 2k dxf< < [ W\ n dx)i = (r±- [ \w\ 2k dx)* 





3rl 






1 1 ^r' 1 



|w| 2 *W)s 



Then on the larger domain we have 



1.48) <(i^__/ | w | 2 ^)i = ( J|4 )i( _L / \ w \ 2k dx)i 



\Br\ 


1 


\B r ,\ 


\B r \ 



r\ J B 



which is the domain over which we may apply the Sobolev embedding theorem 
since w has compact support. Now we will deal with the size of the CC balls with 
the following proposition from Nagel, Stein and Wainger, [29J and use 



Proposition 1.5. Given a bounded set U C R n there exist Q > n, Rq > 0, Cb > 0, 
such that for every x G U , R < Rq and < t < 1 

(1.49) \B(x,tR)\>C B t Q \B(x,R)\ 

In our use of the proposition, we have t = ^ and use 

Claris \B(x,r)\ \B(x,r)\ r Q 

\B{x,tr)\ \B(x,jr)\ r 

Now estimate the size of the balls with Cb(^)® ■ From fll.48j) 
(1.51) < C s [C B ^) Q/k }(r^ I \V Q w\ 2 dx) 

I \ n r\ J B r 
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Now substitute the definitions of w — vip, and Vo^ = V^vip + Vq^v 

(1.52) =Cs [C B C-) Q/k }(T^ [ \V v^ + VoM 2 dx). 

r l-D-rl J B r 

We notice that 

|VoW> + V #| 2 = |V v| fy 2 + ^ 2 |Vo^r + 2^((V v)vVo^} 
Now apply the Schwarz inequality 

< \V v\ V + ^ 2 |Vo^| 2 + 2[(|V„i>| W|V ^| 2 ) 1/2 ] 
And Young's inequality 

< |V «| V + ^IVo^I 2 + [(|V «|V 8 ) + (M 2 )|V V| 2 ] 

= 2|V t;|V + 2^|Vo^| 2 
Since we are on B(x,r), ip < 1 and |Vo^| < Chg/{t — r ') 2 

= 2\V v\ 2 + 2C HG v 2 /(r-r') 2 . 

We replace this back into (II. 52ft : 

(1.53) < C s [C B (^/ k ](j^ [ [2\V v\ 2 + 2C HG v 2 /(r - r'f\dx) 

r | -Or | J B r 

= 2 Cs[C B ^) Qlk ]{^ J b | V v\ 2 + (C HG v 2 /(r - r') 2 )dx ) 

This allows us to say 
(1.54) 

(nrr/ <2C s [c B (- f n-^ [ \v v\ 2 + {c HG v 2 /{r- r ') 2 )dx). 

\ B A Jb t , r' \B r \ J Br 

Now we will use this estimate from the Sobolev embedding theorem to prove 
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Theorem 1.6. Let u G S 1,2 (R) be a positive, weak subsolution to W.l\) in R = 
{B(x,r) x (0,r)}. Also define 

R' := {-BforO x (0,r')} 

/or < r' < r, < r' < r. m <n, are smooth Hormander vector fields 

and the matrix of coefficients aij(x,t) are L°° bounded, measurable and uniformly 
elliptic satisfying 

for every (x, t) G R and for every £ G M" . Then there exists a constant 7 depending 
on X, n, and R such that 

(1.55) maxii < 7 (-— // u p dxdt) 1/p 



R 



for p > 1 . 



Proof: In order to prove the theorem, we need a lemma modeled after Moser 
[25] . Now for v G S X ' 2 (R), define 

(1.56) H, A v): =w ±-^jj^ rl r,M 

(1.57) M rr (v) := 7-— ■ max / v 2 dx 

\B r \ te[o,T]J Br 

(1.58) D r>T (v) := 1 / / |V<H 2 etedt 

I -B P x (0,r)| JJ^ 

Lemma 1.7. Let u G S^ 2 (R). Let k = 1 + | /or n G N, a = /3 = f T/ien, we 
/iawe for some constants C\, C2 > 0: 

(1.59) H r , iT (v k ) < d M r>T (v)n [C 2 H r , >T (v) + D r ,, T (v)] 
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Proof: 

Notice that 



1 1 n-2 2 „ 
- + - = + - = 1. 

a p n n 



In particular, a and (3 are Holder conjugates. First consider 



v 2k dx = / v 2{ - l+ ^dx = I v 2 v«dx 



And by Holder's inequality we have 



<(/ v 2a dx)«( v -dx)?=( v 2a dx)«(l v 2 dx)- 
'Bj Jb„, Jb„, 



<(/ v 2a dx)«M(v)n\B r \n 

JB r , 

Since the max over all t is greater than or equal to the value for any t. We will 
use Sobolev's embedding theorem on the RHS. We notice now the step mentioned 
earlier. The integral 

(1.60) (f v 2a dx)^ 

JB T i 

is taken over the ball, B r > and we will use the results from above. Use the Sobolev 
embedding theorem on the (jl.60p term in the form of the inequality (jl.54p to 
achieve: 

/ v 2k dx <2C s [C B (^) Q / k ]([ \V v\ 2 + (C HG v 2 /(r-r') 2 )dx)M(v^\B r \^ 

JB r , r JB r 

Integrate in t over [0',cr] to get: 

\B r , x (0,7 J )\H r , tTl (v k ) 
<2C s [C B ^f/ k ]\B r x(0,r')\ 
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x [D ry (v) + (l/r 2 )H ry (v) ]M rtT (v)$\B r \$. 
Now we combine constants and divide the size of the \B r > x (0, r')\: 



H,M) < 2 Cs[CB (l)«A] 



X 



[iVM + (l/r 2 )^(t;)]M r , T (t;)n|S r |-. 

Before we finish with the creation of the iteration inequality we will deal with the 
size of the CC balls with the proposition from Nagel, Stein and Wainger, |29| and 
use 

n m\ \B(x,r)\ \B(x,r)\ r Q 

\B(x,tr)\ \B(x, r v r)\ r' 

Now estimate the size of the balls with Cb{^j)^ ■ 

H r ,, T ,(v k ) < 2 C s [C 2 B (-f+Q/ k ][D rtT ,(v) + (l/r 2 )H r:T ,(v)]M riT (v)Z\B r \%. 



(1.62) H r , y (v k ) < 2 C s [Cl{^f + Q' k \\B r fi[D r>T ,{v) + (\)H r y(v)]M r ^ 



(v) n . 



This proves the lemma. In terms of H r ^ T (v), D r>T (v) and M r ^ T (v) we calculate our 
local bounds (II .2^ and (jl.3p as: 

(1.63) M rT (v) = -max / v 2 dx 

\B r ,\te(o,T)J Br 



(1.64) <s\(-^- + ^—)-^- I I v 2 dxdt 



(r-r') 2 t-t''\B. 



r\ J J R 



(1.65) 
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Then for some C3 depending on A, R and n we have 
(1.66) M rtT (v)<C 3 H r>T (v) 
And 



where we need to replace the r' with a r 



1 l(0,r)| ox2/ C\ G , w/ , 



D rT( (u) < — ; ry - 8X 2 (- — + ) // v 2 dxdt 

• K J - \B r x (0,r')| |(0,r)| V (r-r') 2 r-r' ;,/ 



" |S r x(0,r)| |(0,r')| l (r-r') 2 t-t' , JJ r 

t (r — ry t — t 
gives for some C 4 dependent on A, n and R 

(1.67) D r , y (v) < C 4 H r . T (v) 

So substituting these into the results of the lemma f)2.4ip above, we get 



H r , y (v k ) < C x CzH rjT {v)^[C A H r>T {v) + C 2 H rjT (v)} 



dC^C^ C 2 }H r , T {v)-H r , T {v) 



(1.68) =lH r ^v) 
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where 

7 = 2C, s [Ci(^)«+«'*]|B,| 2 '" 

This shows that for a positive, weak subsolution v in R that is square integrable, 
then v k is square integrable in R! . We can then repeat this exponentiation by k to 
achieve the max of v in a smaller domain. For example, take 



?W2 



^1.69) v v := f 



a positive, weak subsolution in R where p v = pok u , p v > p > p' > 1 and k = 1+2/n, 
n E N is as above. Notice that k > 1 for all n and since p is greater than 1, p v > 1 
for all v. The construction allows — > oo as z/ — > oo. Likewise, when p < as 
will be the case for minu(x,t), then p„ — > — oo. 

Now let us set up the domains for each of the iterations. Each iteration will be 
defined in its own rectangle. 
Let 

(1 ,0) ,.:=4±£ 

which converges to r' as v — > oo. Also let 

1 + r'v 

(1.71) r^^r— 

1 + rv 

which converges to r' and v — > oo. Let 

(1.72) = x (0,7V,). 
Notice that since z/ starts at 0, we have 

(1.73) r = r TT -^ = r 
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and 



1 + t' 

(1.74) n = r < r 

1 + T 



which gives us smaller r v as v — > oo. Also notice that 



f ,^r^ 1 + r -0 

(1.75) ro = r ___ = r > ri . 



Likewise, 

(1.76) R = B x (0, r ) = r) x (0, r) 



and R u — )■ /?' as v — > oo. 
Set 

(1.77) H v := H rvTv [v v ) = t^i 1 1 v^dxdt. 

\ K v\ J J R v 

Now we will show the iteration. Take for example, for any positive value of u 

(1.78) H u+l = H rv+UTv+1 {v u+l ) = [ [ v^dxdt < 7 ff* 

l-nV+il J JRv+i 

To arrive at a constant for 7 we need to go back to the local bounds calculated 
earlier for higher powers of p. Using the same conventions for H u = H TvTv we have 

M v+ i = — -max I v u+i dx = — --^max I v Pok " +1 ^ 2 dx 

\B v +i\ t J Bv+1 \B v +i\ \R V \ * Jb v+1 



(r„-r v+l y t v -t u+1 \B v+l \ 



(1.79) < 4p, +1 A( - f + —-—){C B {^) Q )r v H v 
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Where Chg was the constant from the bounds on the horizontal gradient. Cb was 
the constant from the bounds of the CC balls. 



Du+i = | p 1 ■ 1 1 \V v v+l \ 2 dxdt = 1 j^j ! I \V v u+l \ 2 dxdt 



v + l 



(1.80) < A\\ Pv+l f/{2p u - 1)( C2rg + ){C B {^Jf)H v 

(r„ - r v+1 ) 2 r v - r v+1 r u+1 

Now substitute (jl.79p and (jl.80p into the lemma above. This gives 



< [AX\ Pu+1 f/{2 Pu - 1)( °\ G + -±—){C B {^)H u + H u ]x 



< [4A>„ +1 )V(2p„ - 1)( °» c + _L_)(C B (J^)«) + l]x 



[W fa C f + )(C(^-)«)r.]ifl* 

Which is 

(1.81) H v+1 < 

with 7 depending on R, n, A only. Now we will use this fact to achieve a maximum. 
First we will iterate down to Ho. Take any H v and notice that 



for some 72 > 7f . So 



H v < 72^-1 



Now use (H u _i) k to iterate again. 



_ M+(u-l)k rrk 2 
— 72 n u-2 



for some 



So we have 



< 72 172 H u-3) 



(^)+(^-l)fc+(i>-2)fc 2 ^fc 3 
72 n v-3 



/ M+(j'-l)fc+(i>-2)fc 2 +...+(l)fc 1 '- 1 „-fc" 
— 72 n 



< 73 ^0 



73 > 72 



tt ^ „,fc" rrfc" 
#1/ < 7 3 ^0 



which in terms of the integrals says that 



50 

Then taking k u roots: 

1 f f , „./n,i.K . ..vl/fc" 1 



(1.82) [rBl \v Po/2 f dxdt) 1 <^— v po/2 dxdt. 

We now use the following fact from L p theory: 



lim ( / v p dxdt) 1 ^ = \ \v ||oo = ess sup w. 
p^°° Jr r 



So on the LHS of (11.82]) we have 



maxv = lim(-^- T // \v Po/2 \ k " dxdt) 1/k " 



Rl u ^°° \Ru\ JJR V 

So we conclude: 

(1.83) maxM < J J u p dxdt) 1/p . 

We will consider the case for p < in order to achieve a similar result for 
the minimum of t> p °/ 2 in R', where v is a weak supersolution. So assume that 
Pu < Po < where — p is the same value as above. Let p v = p k u /2. Then 
limpj, = — oo. Notice also that at v — 0, we have p u= o = —Po/2 as before. Each of 

v— >oo 

the calculations are the same as the < Pq < p v case except that we have 

H u+l > jH v 

since v is a supersolution we have: 

/ / 4>v t dxdt + / / Xi(f)(ciij(x, t)Xjv) dxdt > 0. 
J J j J=1 </ J 

Then by a similar argument we have 

(1.84) H v > ^H k " 
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and taking k v roots the RHS of (jl.84p tends to the minimum of v Po/2 in R', 
limfr^-r It \v Po/2 \ kv dxdt) x/ku = mmv 



JJr u r> 



where p < 0. 
So we conclude: 



",■(7777 / / " ''(Ivdt) 1 ^ P < minii 



Now it is only left to establish an inequality that will bridge the gap between 
positive and negative po, 



v pokV/2 dxdt < 7— |— // v~ mkV/2 dxdt 



R-\ \R ■ i j jr 



1.5 Bridge step 

Above we established appropriate estimates on the maximum and minimum 
of a positive, weak solution to (10. ip in a rectangle R. We turn to the gap for 
— p < < p . Our earlier results gave us (I1.83P 

maxu p < 7 / / vPdxdt 

RI JJr 

so we start with v = u~ e which was proven to be a positive subsolution under the 
p < case. Applying the results of (I1.55P above we get 



maxn e < j(J J u €p dxdt) 
Now taking the — 1/e root and setting p = 2 we find 



mmu > 7~ 1/e ( J j u~ 2e dxdt)~ 1/2e 
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Let 7~ 1 / e = 7 X and use e for 2e 

(1.85) mmu > ji{ J J u~ e dxdt)- 1/e 
Introduce the notation for some domain D C R 

(1.86) M(p, £>) := (r^-r / u p dxdt) 1/p 
and use this for (II. 85ft which gives us 

(1.87) M(-oo,R') = rmnti > 7i(^ u- £ dxdt)~ 1/e = ^M{-e,R). 
We also want to establish for some constant 7* 



1.88) max u = M(oo, iT) < 7*( // u e dxdt) 1/e = 7*M(e, 



So we let R* C R where the boundary of i?* does not touch the top of R . Now 
apply the results of (I1.55P to u in R* with p' > 1, then for some 7': 

(1.89) maxii < ^'M{p',R~) 

R* 

Recall that we proved the local bounds: 

|V f \ 2 dxdt <C 1 J J v 2 dxdt 
and 

max / v 2 dx < Co v 2 dxdt 
te(o,T)7 J J 

for < p < 1, in the form v = u p ^ 2 (—t,x). We create an iteration argument like 
the one above. This gives us, for some 7 the inequality 

(r^r J I U pk dxdt) 1/Pk < 7 (^ J J U P dxdtf lp 

like before. Now let pk = p' and p = j-. This gives 
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1 



u p 'dxdt) 1/p ' < 7 (^ T ff u p ' /k dxdt) k/p ' 



\ R I JJr- \ R \ J JR 

The issue here is that the p' will stay fixed and we will iterate the RHS to achieve 
£7 — > e with larger Using the k = 1 + 2/n from the previous iteration, we set 
p' = *±1 and = ^^r-. This § ives 

M(p',R-)< lv M(p V} R) 
for some 7„. Now let ^ be the smallest integer such that 

(1-90) = l kJ ^ < ' 

We find 

(1.91) M(p', R-) < 7 M(e, R) 

Combining (jl.9ip and ( I1.89P gives 

max u < 7'M(p', iT) < 77'M(e, i2) 
which proves ( 11.88P with 7* = 77'. Take this along with 



M(—e,R) < 7minw 

proved earlier. To finish the proof of the Harnack inequality, we need to bridge the 
gap between M(e, R*) and M(— e, R'). In order to do this we will use the results of 
Hugo Aimar [T]. 



1.6 Lag mapping 



First, some definitions. 
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Definition 1.8 (Lag Mapping). Let R(x ,r ) be a fixed ball in the space of homo- 
geneous type (X, d, /j). We use X = IR n+1 . We shall say that a function 



is a lag mapping if there exist three constants K it i = 1, 2, 3 such that the inequalities 



hold for every x G R(x , r ) and r G (0, r ]. 

In the following we will consider rectangles 

R((x, t),r) = R(x, r) = B(x, r) x (t — r, t + r). 

In the definition of BMO we shall use non-negative functions h of a real variable t, 
which are continuous, non-negative, h(t) — if t < 0, h(t) is increasing for t > 0, 
h(t + s) < h(t) + h(s) and e~ eh ^ is an integrable function on (0, oo) for every e > 0. 

Definition 1.9 (BMO condition). Let f be a real-valued measurable function on the 
space of homogeneous type (X, d, /i). We shall say that f satisfies a BMO condition 
with lag mapping T on the rectangle R(xo,ro), with respect to h(s) = y/s if there 
exists a real function C(x,r) on R(x ro) x (0, Rq] such that 



T : R(x , r ) x (0, r ] X x R+; T(x, r) = (£, p), 



(1.92) 



d cc (x, < Kir, 



(1.93) 



K 2 p <r< K 3 p 



(1.94) 




(1.95) 
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for some constant N(f) and every (x, r) G R(xq, r ) x (0, Rq] . Let BMO(xo, r , T, h) 
be the class of all such functions. 

The main result of the paper and the one we will make use of is the 

Theorem 1.10 (John Nirenberg for BMO with Time Lag). Let T be a one to one 

lag mapping on R(xo,ro) with the following property: there exists < 7 < 1 such 
that for every r < r$ 

(1.96) R(x ,-fr) x (0, 7 r] C T(R(x ) x (0,r]). 

Then there exist two constants 5 and C depending on Ki, A, 7 and a lag mapping 
S on B(xq, 5ro) (S(x, r) = ((, r)) such that if f G BMO(xq, ro, T, h) and u = e~f , 
the inequality 

(1-97) (m R ( x ,r)( u ~ € )) ■ ( m R«,r)( u€ )) < C 

holds for some e > 0, depending only on N(f), and every (x,r) G R(xo,5r ) x 
(0, 5ro). Moreover, the mapping S can be given explicitly in terms ofT: 

r r 
S(x,r) = (f(x, -);r)p(x, -)), 
77 77 

where £ and p are the components ofT and r) is a constant depending only on Ki, 
i = 0,1,2,3. 

We will use this in the following way. The bridge part of the theorem is from 
(I1.97P and takes the form 

(1.98) -—. 1 . ., / / u~ e dxdt ■ 1 1 u'dxdt < C 

\R({x,t),r)\JJ R{{xAr) \S(R)\JJ S{R) 

Where 

S((x, t), r) = (x, t - 2r 2 /7] 2 ),r) 
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Now we want to show that this theorem applies to our setting. First our space 
(X, d, fi) is a space of homogenous type with 

X = R rt+1 

d((x, t), (y, s)) = d cc (x, y) + \t- s\ 1/2 

and fi = (n+ 1) -dimensional Lebesgue measure. We will use two notations: (x, t) G 
X C R n+1 or x G X C R n+1 depending on whether we want to draw attention to 
the time dimension of a point. 

The constants K^p < r < K%p are both 1. Now we need to show that the 
doubling property holds. We want an A such that 

< fi{R{x,2r)) < A/i(R(x,r)) 
holds for r > and every (x,t) G X. An earlier proposition |1.5j gives us 

\B(x,tr)\ > Ct Q \B{ x,r)\ 
for every x G U bounded, r < R and < t < 1. Sosett = l/2 and use 

\B(x,K)\>C B (±f\B(x,r)\ 
\B(x,-r)\ x 2r > C B (^) Q \B(x,r)\ x 2r 
\B{x,\r) x (^, T -)\ > C B (±f +1 \B(x,r) x (-r,r)\ 

\R{xJ-)\>C B {\)^\R{x,r)\ 

Which gives the desired inequality for A > ^■(2) < ^ +1 

We will use the lag map T:K n xf x[0, 00) — > M. n x IR x [0, 00). This map acts on 
pairs from M™ +1 and R + in the form T((x, t),r) where (x, t) is a point in space-time 
and r is a radius for the rectangle in IR n+1 . Define T((x, t); r) := ((x, t — 2r 2 ); r) = 
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((£, s),p). We want to check that T is a lag map. Notice that 

t), (£, p)) = rf cc (x, x) + |t - pi 1 / 2 = + |t - (t - 2r 2 )| 1 / 2 = rv^. 
Last we check that there exists < 7 < 1 such that for every r < R 

(1.99) R(x ,jr) x (0,tt] C T(R(x ,r) x (0,r]). 

i?(x , 77-) x (0, 77"] = B cc (x, yr) x (t — 77*, t + 7r) x (0, 77"]. 

T(R(x , r) x (0, r]) = 5 cc (x, r) x (t - 2r 2 - r, t + r - 2r 2 ) x (0, tt] 

Clearly, B cc (x^r) C B cc (x,r) and (0,77"] C (0,r] for every < 7 < 1. And 
(t — jr, t + 77") C (t — 2r 2 — r, t + r — 2r 2 ) for 7 < 1. 

1.7 Estimates on log(u) 

Now we calculate estimates on v := — logw. We will use the results of these 
estimates to show that v G BMO(xq, Tq, T, h) in order to use the results ( 11. 97ft 
of Aimar which will conclude the proof. We will consider positive supersolutions 
u > e > in B{x,r) x (— r < t < r). Let v(x,t) := —log(u). Recall from the 
definition of weak supersolution to 

m 

u t = - 2J X*(a ij (x,t)X j u) 
in the region R defined before: 

(1.100) J J (j)u t dxdt + ^2 J J (a ij (x,t)X j u)X i (j) dxdt > 0. 

We will use = ip 2 u~ l where (f)(x,t) > is of compact support in the x variable 
and ijj is independent of t and of compact support. 
Then we have with 

v t = {-log(u)) t = -— 
u 
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XiV = Xi(-log(u)) = 



u 

-1 „.-2v „,i,2 



= I^X^u' 1 - u- z XiUil) 



(1.101) / 1 r 2 -^ • < 2rV r// ; .r/V // > r/.rr// 

< ^^V 2 ,aV M > dxdt > 0. 
And substituting for u 4 and VqV 



■{ip 2 v t )- < V (ip 2 ),aV v > - < ip 2 V v,aV v > dxdt > 



R 



R 



ip v t + < V (ip 2 ),aV v > + < ip z V v,aV v > dxdt < 
Now consider the integral in t from t\ to t 2 in the interval: 

' B(x,r) Jti J B(x,r) 

And taking to the RHS the second integrand 

-t 2 



[ ip 2 v dx\\\ + [ [ 2^ < Vo-0, aV w > +ip 2 < V v, aV v > dxdt < 0. 

J B(x,r) JU J B(x,r) 



^ 2 vdx\ l t \+ I [ ip 2 < V v, aV v > dxdt 

x,r) J t\ J B(x,r) 

< -2 / [ < V ip, aV v > dxdt. 

Jt\ J B(x,r) 



(1.102) f ip 2 v dx\ll + [ [ ip 2 <V v,aV v > dxdt 

JB(x,r) Jt x JB(x,r) 



< \2 f f ip < V ^,aV v > \dxdt. 

Jti JB(x,r) 
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Now we look at the second term in the integral and use Schwarz' and Young's 
inequality: 

2 / / ip < Vo^, aVof > dxdt 

Jt! JB(x,r) 

< 2( [ [ ip 2 < V v, aV v > dxdt 

Jt! JB(x,r) 



X 



[ [ <V v,aV v> dxdt) 112 

Jt x JB(x,r) 



(1.103) < 2{\ [ [ ip 2 <V v } aV v > dxdt 

4 Jt! JB(x,r) 

+ / / < Vof, aVo^ > dxdt) 

Jt! Jb{x,t) 

When we combine (11.1021) with (11.1031) we get: 

(1.104) [ ip 2 vdx\ t t 2 1 + - [ [ ip 2 <V v,aV v > dxdt 

J B(x,r) 2 J tl Jb(x,t) 

< 2 [ [ < V ip, aV ip > dxdt 

Jt! J B{x,r) 

Now we make use of Jerison's Poincare inequality for Hormander vector fields 
We will point out that Jerison assumes v is C°°(B(x, r)). In order to improve this 
result for functions v G S 1,2 (B(x, r)), we will use the fact that the smooth functions 
are dense in S l,2 (B(x, r)). I.e. given a v G S 1,2 (B(x, r)), there exists a sequence of 
smooth functions on the ball, {v n }™ =l G C 00 {B{x^ r)), such that 

lim v n = v 

n— >oo 

and 

lim Vof n = Vo^ 
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in the S 1 ' 2 norm. We will let Q be a neighborhood of the closure of the unit ball 
B(x, 1), in R n . 



Lemma 1.11 (Poincare inequality). Let {X i }™ =1 , m < n, be smooth vector fields 
satisfying Hormander's condition. Then there is a positive constant C p and a radius 
ro > such that for every x G B(x,l) and every < r < ro for which B(x, 2r) C Q 

(1.105) / \v -V B \ 2 dx <C P r 2 \V v\ 2 dx 

JB(x,r) JB(x,r) 



holds for every function, v G S 1,2 (B(x,r) x (0,r)) with fixed t. Here Vpit) = 
fv(x,t)dx. 

Now we construct a test function whose existence is guaranteed by Lemma (II. 2ft 
|10j . Let ijj(x) be defined as follows: 



ip(x) 



1 if dcc(x, 0) < r 
if d cc {x, 0) > 2r 



dx \t + ^l I \y v\ 2 ip 2 dxdt < 2X(t 2 -h) [ \W^\ 2 dx. 

J tl J B(x,r) JB{x,r) 



We note that (Vo^l < Cnc/r 2 - Now we recall the bounds on aij(x,t), integrating 
in t on the RHS since ip is independent of t and use (11 . 104j) . 

^dx \H + ^ I 2 

I B(x,r) 

We observe that we defined ip(x) = 1 on B{x,r) so this leads to 

>B(x,r) ^ Jti JB(x,r) J B(x,r) 

Substitute the value Vb = n^n, 31 f B r x r \ vdx and divide by (ti — ti) 

■t 2 



1.106) f vdx\il + ^-[ [ \V v\ 2 dxdt <2X(t 2 -t l ) [ \V tp\ 2 dx. 

JB(x,r) 2A J t JB(x,r) J B(x.r) 



v B (t 2 ) yi) )|B(g>r)|+ 1 r /■ | VoV |2 M < 2A /" | Vo ^|2^ 

l J 2 — <^ Jt! JB(x,r) JB(x,r) 
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Now use the Poincare lemma 

{ v B (^-v B (h) mx ^ + c^_ r<> r (v _ VB?dxdt £ 2A f 

V>2 — tl) lA{t 2 — ti) J tl JB(x,r) J B(x,r) 

On the RHS we have |VoV'| < 1/V so we get 



(*2 - ^l) 2A(t 2 - £l) 7 tl 7B(x,r) ' r2 JB{x,r) 

Now divide through by \B(x,r)\ 

VBi t - If' 1 ' + 2|g(x, S h - tl) C Lj v - VB?dxdt - wBw^ r)] 

Take t\ — >■ We note that Vb need not be differentiable although we will assume 
so for the present argument. We present a remark (11.121) at the end of the proof over 
an alternative technique using the Steklov average. This allows us not to assume 
the t derivative of Vb exists. 

(1.107) ^ + ° pT \ - [ (v(x,t) - V B fdx <^=:m 



dt 2\B(x,r)\X J B (x,r) ' ' r2 

Let w(x, t) = v(x, t) — mt and W = Vs{t) — mt. Then 

dW dV B 

— — = — m 

dt dt 

and 

v(x, t) -V B = w(x, t)+mt — (W + mt) = w(x, t) - W 
Substituting this back into (I1.107P we get 
dW C r 2 f 
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Let 

B s (t) = {x e B(x,r) : w(x,t) > s} 
Then (s - Wf < (w - W) 2 , hence 

<W C P r2 f / r,rN2 , 

1.110 — < — , , — / (s-MOtfe 



We integrate from ti to t 2 : 

/■*2 AYA7 n r 2 

(1-113) / -tt(s - W)- 2 <- -. ^ P - \B s (t) x (t lt t 2 )\ 



ti 



dt y ' ~ 2\B(x,r)\X l 



1 1 x C P r 2 . 

< - 1|D , m X |^(*) x (h,t 2 )\ 



(s-W(t 2 )) (s-Wih))'- 2\B(x,r)\\> 



(s-W(t 2 )) (s-Wfa)) - 2\B(x,r)\X 
We defined W(t) = V B (t) — mt and from (I1.107j) we see that 

(1.114) < 2 4 = m. 

y ' dt ~ r 2 

So we can say W(t) < on (tx,t 2 ). So that leads us to say: 



11 1 1 

> — - > 



s ~ (s - W{t 2 )) ~ (s - W(t 2 )) (s - W{t x )) ' 
Hence 

2\B{x,r)\\ c 3 



(1.115) \B a (t) x (t u t 2 )\ < 



Cpr 2 s r 2 s 
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Set 

R((x,t),r) = B(x,r) x (t-r 2 ,t + r 2 ) 

So we can say that 

(1.116) n{(x, t) G R((x , t ),r):w> s} < -|- 

f s 

for every s > 3mr 2 . We let /i be the n + 1 dimensional Lebesgue measure on the 
rectangle. We will prove this for w since we defined w = v — mt. So 

(1.117) / vfej - V B ) + dfi(x,t) < Vw+dfi(x,t) 

jR((x ,t ),r) J R((x ,t ),r) 

Use the Lebesgue-Stieltjes integral: 

I V w + dfi(x, t) = — / s~ l l 2 dii(s) 

J R((x ,t ),r) 2 J 3mr 2 

And integration by parts 



= -\Us)s^ 2 \f mr2 - \ S^ 2 ^ S )ds] 
A 1 J3mr 2 

with the boundary terms equal to 0. 

oo 

s~ 3/2 n(s)ds 

J 3mr 2 

Now we use the bounds on jj,: 



-y poo 

4 



4 J 3mr 2 r 2 s 



;i - n8) -l? 2 L^ 2 



oo ^ 

ds 
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Which converges since 

C I —,ds < C / — ^ds < oo. 



3mr 2 s5//2 JR s5 ^ 2 

This leads us to say 

/ y/(v(x,t)-V B ) + dfl(x,t) < C 4 

JR{(x,t),r) 

Now notice that for any slice of time t\ to we have \t\ —t%\ = 2r 2 . If we replace 
t with 2t\ — t in the above we find 

h < 2*i - t < t 2 

-ti <-t<t 2 - 2*i 
h>t>2t 1 -t 2 = 2*i - (* + r 2 ) 

and that gives 

|(*i) - [2*i - (* + r 2 )]| = |*o - *i - r 2 | = 2r 2 

This is the correct size of the time slice with *i as the upper and t 2 as the lower 
bounds. Now we consider 2*i — * G (*2,*i)- We take the argument over at ( 11.1 14ft . 

dV B 2A _ 
dt ~ r 2 

and we integrate from t\ to * to get 

(1.119) Vb(*i) -V b < m[{2h - *) - * x ] < 2mr 2 . 

Now for s > 3mr 2 set 



= {xe B(x, r):V B - v(x, *) > s} 
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which implies V B (ti) — s > v(x, t) Putting these requirements together with fl 1 . 1 1 9 [) 
gives us: 

V B — v(x, t)>V B - (V B (ti) -s)> s- (Vs(ti) -V B )> s- 2mr 2 > 

Set 

R((x, t),r) = B(x, r) x ((2t x - t) - r 2 , (2t x -t)+ r 2 ) 
We use the above inequality with ( 11.1071) to achieve following the integration above: 

fi({x G R : V B -v(x,t)}) < ^. 

This converges by the argument above. This gives us 

f ^/{V B -v{x,t))+d^{x, t) < c 5 

J B(x ,t -2r 2 )x(-r,r) 

where the center of the new ball is t — 2r 2 as expected. 

This shows v = —log(u) belongs to BMO with h(s) = a/s+ and lag mapping T. 
Aimar's main theorem then gives the desired result 



( f u- £ dfx) ■ ( f u £ dfx) < Cfi(R) 

jR((x,t);r) J R((x,r);r) 

where r = t — 2(r/rj) 2 . In the above 



77 = [K {1 + 2[Kl(2(K^ + 2K 1 )K 3 Kz 1 + 5K 2 K^ + K x ) + l]#o)] _1 
This bridges the gap if we choose from (jl.90p p v = e. 

Remark 1.12. We begin our discussion of an alternate proof with recalling equation 
[LIM) for a weak solution u E S^R), R C W l+1 . 

(1.120) J J cj)u t dxdt+Y^ J J (a ij (x,t)X j u)X i (j) dxdt > 0. 
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From M.lOfy) we proved the inequality 

(1.121) ! ^J{v{x,t)-V B )+ dfJL(x,t) < c 4 

JR((x,t),r) 

Now we will alter the argument from the beginning with the Steklov average. 
Definition 1.13. Steklov average 

For a positive, weak solution u G S 1,2 {R) to ( fO. 1\) we define 

u h = - J u(x,T)dr 

Observe that if u is a weak solution to W.l\) then so is Uh- Also, we have that uy t is 
differentiable in t so we can use the following form of the weak solution: For each 
t G [£1,^2] C [0, T] we have 

(1.122) J J (j> d t Uh{x,t) dxdt + J J {aij(x,t)XjUh(x,t))Xi<f> dxdt = 

for every <p(x, t) G C°°(R) for each fixed t. 

We also have the following results from L p theory concerning the Steklov average 
[23]. 

i. Ifu G L p (n x (0,T)) then u h ^ u in L p (Vl x (0,T - 5)) as h -> 0+ for 
< 5 < T. 

ii. IfVu G L P (VL x (0,T)) then Vu h -»■ V« m L P (Q x (0,T - 5)) as h + 
for < 5 < T. 

Here Vu denotes the Euclidean spatial gradient of u. Notice that the analogous 
result for the horizontal gradient of u follows from (ii). 

We may continue from ^1.122) to achieve the same bounds on log(u). Then, at 
the end of the proof of the bounds on log(u) before the application of the modified 
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John-Nirenberg theorem of Aimar we will have from U.98\) 

(1.123) , „ 1 r / / u- £ dxdt ■ — - / / uldxdt < C 

1 ; \R((x,t),r)\JJ mX;tM h \S(R)\JJ S{R) h 

Now we may apply the results from the LP theory of Steklov averages: 

lim —— r — -- / / u7 e dxdt ■ . - . — , / / uldxdt 

h->o + \R({x,t),r)\JJ R{{Xit)ir) h \S(R)\JJ S{R) h 

^ ^ u~ e dxdt ■ / / u e dxdt < C. 



\R((x,t),r)\ JJ R{{xA ,r) ' \S(R)\jJ S (R) 

1.8 Holder Continuity 

We now show that weak solutions in a bounded domain are Holder continuous. 
This section follows closely the approach by Moser [25J. Let us construct appro- 
priate domains to examine the continuity. First let u be a positive, weak solution 
to 

in 

u t = - ^ X*(a ij (x,t)X j u) 

in 

R = B(x,r) x (0,r) 

and sub domains 

R+ = S(x,r/4) x (7r/8,r) 
iT = B(x,r/A) x (r/8,r/6) 

We will use the distance t), (0, 0)) = max{d cc (a:, 0), yjtj} for x G R™ and 
i Gl, Also, we will form ann + 1 dimensional vector x' = (xi,X2, ...t) G R C M n+1 . 
With this extra notation we will now restate the distance d((x, t), (y, s)) as 

d(x',y') := d((x,t), (y,s)). 

Then we have for a positive, weak solution in R, u(x,t) := u(x') a corollary to the 
main theorem: 
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Corollary 1.14. (Holder Continuity of weak solutions) Let u e S^ C (R) be a weak 
positive solution to 

m 

Ut = -^2 X*(a ij {x,t)X j u) 

in a domain R. Then using the notation from above, there exist constants C > 0, 
and < a < 1 such that for every x' , y' e Rl C R, x' ^ y' , we have 

(1.124) \u(x')-u(y')\<Cd(x',yT 
Set the average of u in R~ : 

:= 1 1 u(x,t) dxdt < maxw 

\R I J Jr- r ~ 

We proved earlier in the main theorem that 

(1.125) minw > 7"V _1 = 7 _1 ^^ / / u(x,t)dxdt 

R+ \R I J Jr- 

For p = 1. For any positive, weak solution u in a bounded domain i?, let 



M = ess max -u 

R 



m = ess mmu 

R 

M + = ess max u 

R+ 

m + = ess min-u 

R+ 

Then M — u is a non-negative weak solution to 

m 

«t = - 5^ Xt(aij(x,t)Xju) 

since it is const — tt > 0. Likewise with u — m. The mean values of M — u, u — m 
in i?~ are 

j^j- ^ (M - u) drrdt = M-n~ 



(>!) 



\R- 



(u — m) dxdt = fj, — u 



r- 



Apply (11.1251) to both 



(1.126) min (M — u) = (M - M + ) > 7 -1 (M - //"). 
Since (M — u) is minimized when u is maximized in R + . 

(1.127) min (u — m) = (m + — m) > 7~ 1 ( / u~ — m). 

Since (u — m) is minimized when u is minimized. Let the oscillations u = M — m 
and u> + = M + — m + represent the difference between the max and min on a given 
domain. Adding fll.!2f)|l and fl!.127j) we find that 

(1.128) (M - M + ) + (m + - m) > 7~ 1 ( / u~ - m + M - 
implies 

(M - m) - (M + - m + ) > 7~ 1 (M - m) 

u — 7 _1 w > uj + 
w(l -7 _1 ) > w + 

So we can say 

'Y — 1 

(1.129) w + < (- )w = du 

7 

for some < 9 < 1. Let x' £ R' C R with a fc-neighborhood also contained in R '. 
So for any x', y' G R' such that d(x',y') < k then the oscillation of u in i?' is w. 
Also we have if d(x', y') < | then the oscillation of u in R! is 9uj. Finally this allows 
us to say if d(x', y') < then the oscillation of u in R' is 9 u u. 
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Let R' C R be any region. Let x' G R' with a ^-neighborhood contained wholly in R. 
Then for x r , y' G i?' satisfying d(x', y') < k, choose v such that -^fi < d(x', y') < Jr- 
Then we have Jj < ^d(x',y'). Putting these together we have 

(1.130) \u(x') - u(y')\ < e»co < = < [~d{x\y'Wu 

= [(^-) a oj}d(x',yT = Cd(x',y'r 

with 

_ -log{0) -logj 1 ^) 
a ~ log(A) ~ log{A) 

and < a < 1. This proves the Holder continuity with respect to the CC metric and 
so our second main result. We notice that with the above result and a proposition 
of Nagel, Stein and Wainger [29j we have that the solutions to (10. ip are Holder 
continuous with respect to the Euclidean metric. They proved that there exists a 
constant c > and a constant e depending on the number of commutators it takes 
to span the space such that 

c\x -y\< d cc (x,y) < c^\x - y\ e . 

So with HI . 1301) we have 

\u(x,t)-u(y,s)\ < Cd(x',y') a < c' l C\x-y\ ea . 
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CHAPTER 2 
Nonlinear Generalization 
2.1. Introduction 

Now we will consider a more general non-linear case of which the above results are 
a special case. We will examine the operator 

m 

(2.1) Lu = Ut + ^2x*A(x,V u) 

3=1 

And corresponding solutions to 

(2.2) Lu = 

Where A : M. m — > M. m is a measurable function, and X* is the formal adjoint of Xj. 
Let u G S^^R™ 4 " 1 ). We want to consider weak solutions u to (12. 2p in a rectangle 
R((x,t),r) = B cc (x,r) x (0,r). Here, r G IR, x G M n . By B cc (x,r) we mean a 
Carnot-Caratheodory ball B(x,r) := {y G M. n : d cc (x,y) < r} where d cc is the 
natural control metric associated with the vector fields {Xi}™ =1 . In the following 
we will refer to d cc (x, y) as d(x, y). By a weak solution to Lu = in R((x, t),r) we 
mean 

(2.3) / 0«t <ir<it + / ^'(^j Vou)Xj(j) dxdt = 0. 

jR{{x,t),r) j=1 JR((x,t),r) 

for all G Cg°(S(x,r)) with <j>{x) > 0. 

We will consider the following structural conditions on the function A which we 
will denote (S). 



Aj(x, V u) <C\V u\ 
A|V m| 2 < Aj(x, V M)X i M < A|V m| 2 



for some constants C, A, A > 0. 
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Notice that in the parabolic case we considered for the Moser iteration, 

Aj(x,Xju) = a,ij(x,t)XjU. 

We will now develop minimum and maximum inequalities for weak solutions to 
( 12.21) in a rectangle. 

Theorem 2.1. Let R be a rectangle, R((x,t),r) = B cc (x,r) x (0,r) in M. n+1 . Let 

m 

Lu = u t + y~]XjA(x, V w) 
i=i 

and consider positive, weak solutions u G S 1,2 (R) to Lu = in R. Here {Xi}"] =l , 
m < n, are vector fields satisfying Hormander's condition for hypoellipticity . Under 
the assumption of the structural conditions (S) we have the following 

(2.4) max u < 7i(^ / u p dxdt) 1/p 

(2.5) (-- / u p dxdt) 1/p < 72 min u 

\R\Jr ri 

for some constants 7i, 72 > depending on R, n, A ; A, and C . 

With the proof of the preceding two inequalities, we will move to the corre- 
sponding bridge step 




as in the parabolic case above. Together, these three inequalities will give us a 
Harnack inequality for weak solutions to ( 12. 2 j) in R. 

2.2 Non linear Cacciopoli Inequality 

We will begin the proof of a Harnack inequality in the same manner as above 
with the Cacciopoli inequality. For convenience, we will let R = R((x,t),r). 
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Theorem 2.2. Local bounds for weak solutions to Lu = 

With the hypotheses above, including the structural conditions (S), let u G 
S 1,2 (R) be a positive, weak solution to 

(2.7) Lu = 

in R = {B{x,r) x (0,r)}. Also define R' = {B(x,r') x (0,r')} for r > r' > 0, 
t > t' > 0. Then there exists constants C\, Ci depending on \, n, and R such that 

(2.8) / \V Q u\ 2 dxdt < Ci / u 2 dxdt 

JR' Jr 

and 

(2.9) max / u 2 dxdt < C 2 / u 2 dxdt 

te[o,a] J B (x,p>) Jr 

Where the max is the essential supremum over this interval. 



Proof: 

Consider a positive, weak solution u to Lu = 0: 

(2.10) (pu t dxdt + y~] / Aj(ar, X7 u)Xj(j)dxdt = 

Now we let the test function = uip 2 where ip is a function compact in the x 
variable and independent of t. Also, we construct the test function to give ip — 1 
on i?'. Substituting gives: 

/ uvp 2 u t dxdt + / Vom)^(m-?/' 2 ) dxdt = 

Notice that = uu t and Xj(uip 2 ) = XjUip 2 + 2uipXjip. This gives us: 

(2.11) f dxdt + jr f A 3 (x,V Q u)X jU i) 2 dxdt 
Jr 2 =1 ifi 



74 



m „ 

+2V / A j (x,V Q u)wpX j il>dxdt = 
j=1 JR 

(2.12) / — —ip 2 dxdt + Y] / A j(x, VoujXjmp 2 dxdt 

m „ 

= — 2 / Aj(x, Vou)uil)Xjip dxdt 
3=1 Jr 

Now we apply the Schwarz and Young's inequalities to the RHS of ( 12. 12ft . 

( 2\ m « 

— -if) 2 dxdt + 2J. / Aj(x,Vou)XjUip 2 dxdt 



< - / {A^x, Vqw)! 2 ^ 2 dxdt + 2 I u 2 \V ip\ 2 dxdt 



1 

2 JR J 

and recall the bounds on Aj from (S): 

(« 2 )* ,2 



/ ^-^^ 2 dxdt + X [ \V u\ 2 ip 2 cfedt 



C 2 

< — / iVottlV 2 dxdt + 2 I \V i/j\ 2 u 2 dxdt 



R JR 



2 

And adding like terms we arrive at: 

(2.13) / dxdt + (A - — ) / |V «|V tfetft < 2 / |V ^|VtM 
Jr 2 2 Jr Jr 

Notice that |("0 2 M 2 )t = ipiptu 2 + uutip 2 . And use this fact by adding ipi^tu 2 to both 
sides: 

(2.14) [ -(ip 2 u 2 ) t dxdt+(\-—) [ |V «|V dxdt < 2 [ u 2 (|V </f+?M) dxdt 
Jr 2 2 J R J R 

Now integrate this over R a = {B(x,r) xO<t<a<r} and drop the second 
integrand on the LHS. It is > 0. We'll choose a so that 
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(2.15) / u 2 \ t=a dx > -max / u 2 dx. 

Jb(x, p >) 2t£(a,a) JB(x,r') 

This gives us from (12. 14ft : 

(2.16) / -(kV) ( M<2 / u 2 {\V ip\ 2 + \ipipt\)dxdt. 
Jr, 

Now consider the 



R<r ^ J Ra 



(2.17) max / u 2 dx < 2 u 2 \ t=a dx 

te(0,o-) J B ( x ,r') JB(x,r>) 

which we constructed above in fl2.15p . tj) = 1 on R' and t is fixed, so 



max / u 2 dx < 2 v 2 \ t=a dx 

te(0,a) JB(x,r>) J B(x,r') 



< 2 u 2 ip 2 \t =(T dx 

JB{x,r') 



< 2 u 2 ip 2 \ t=(7 dx 

JB(x,r) 

Now recall (I2.16P integrated over R a . We will combine that with the | factor as 
well as the fact that the above is twice the (12 . 16[) inequality. Putting these together 
gives: 

(2.18) max / u 2 dx < 8\ u 2 (\ V ^| 2 + \^t\) dxdt. 

*£(0,<t) J B (x,r') JJR a 

Then use the non-negativity of the RHS and take the integral over R for 

(2.19) ~ 8 // m2 (I Vo ^| 2 + \Wt\)dxdt. 

Now we use as result on the bounds of the horizontal gradient of test functions in 
Citti, Garofalo, and Lanconelli [10J in order to attain a constant independent of if). 
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Lemma 2.3. Existence of Carnot-Caratheodory test functions. There exists r > 
such that given a metric ball B(x,r) CC R, with r < ro and < r' < r, there 
exists a function ip 6 C^°(B(x,r)) such that < ip < 1, ip = 1 in B(x,r') and 
IVq^I — ^~^r- Here, Chg > is a constant independent of r and r' . 



Define ip(x,t) := ipi(t)ip2{x) in the following way, using the cut-off functions of 
Citti, Garofalo and Lanconelli [10J: 



1 if r' < t < r 
if t < 



1 if d cc (x, 0) < r' 
if d cc (x, 0) > r 

and as before we have a bound on the horizontal gradient. This allows us to say 



(2.20) |Vo^| 2 + m < {y^y^^) 
Which gives us 

(2.21) max / u 2 dx < 8( . Crg „ + —^—) II v 2 dxdt. 

This proves (12. 9p . We note that the result is similar to the parabolic case since the 
approach was the same. The difference is that the bounds on the aij(x,t) matrix 
are gone and we find only the bounds from the structural conditions (S). This is 
also true the following proof of the bounds on the derivatives. 
Now we return to (I2.14p 

1 C 2 f F 

-{^ 2 u 2 ) t dxdt + (A ) / |V w| V 2 dxdt < 2 / u 2 (\V ij\ 2 + ijnj) t ) dxdt 

r 2 2 J R J R 
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and drop the first integrand. Integrated over R' gives 
(2.22) (A - — ) / |V «| V dxdt < 2 f m 2 (|V ^| 2 + ipipt) dxdt 



Rt JR 



Again use the fact that if) — 1 on R' 



2 

Which gives us 



C 2 

(A-—) / \V u\ 2 dxdt < 2 / u l (| VoVf + \Wt\) dxdt 

R' JR' 



c 2 

V u\ 2 dxdt < 2(A - — ) _1 / u 2 (| V ^| 2 + |^ t |) dxdt. 



R! 2 j R , 

And using the previous bounds on if), if> t and |Vo^| 

(2.23) / \V u\ 2 dxdt<2(\-— )- 1 ( / ° 2hg xo + —!—) [ u 2 dxdt. 

Jr' 2 (r-r'y t-t' J r 

And the integral on the RHS is taken over R since it is non-negative. This proves 
the result in (12. 8ft . and so Theorem 12.21 is proved. 

2.3 Nonlinear Cacciopoli for u p 

Since we have a Cacciopoli inequality for weak solutions u, now we look at 
powers of u, v := u p . 

Corollary 2.4. The same inequality of Theorem \2.S\ holds, up to a constant, for 
powers of a weak solution u, v = u p . 

Proof: Let (f> = ip 2 u p for p > \ and if) be defined as above. Calculate 

v t = [u p ) t = pvP^Ut 

and 

Xj((f)) = Xj(ip 2 u p ) = 2il)u p Xjil) +pif) 2 u p - 1 X j u. 
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We substitute the definitions of <p and v into (12. 10p : 

m „ 

ifu'fyu^Ut dxdt + V / ^-(ac, V u)(2ijm p X j il; + p^u^Xju) 

,=i J R 



which gives 



m „ 

i) 2 u p {p)u p ~ l u t dxdt + y~] / Vo^ipip^u^Xju) dxdt 

3=1 JR 
m „ 

< -2V / Aj{x^vu)i)u p Xji) dxdt. 
o=i Jr 

We use the Cauchy and Young's inequality on the RHS: 

m „ 

i) 2 u p {$>)u p ~ x u t dxdt + y~] / Aj(x, V Q u)(jnl; 2 u p ~ l Xju) dxdt 

.i=i Jr 



< - I \Aj(x, V u)|V dxdt + 2 [ u 2p \V ip\ 2 dxdt 
2 Jr Jr 

Use the bounds on \Aj\ from (5") and substitute for v: 

1 



i? 2 j B 



ip 2 (v 2 ) t dxdt + A I \V v\ 2 ip 2 dxdt 

1 



< 
~ 2 

Adding like terms to get 



4> 2 {v 2 ) t dxdt + A / |V(HV 2 drrcft 
2 J R 

C 2 



\V Q v\ 2 ij) 2 dxdt + 2 I v 2 \V ip\ 2 dxdt 
Jr 



I C 2 /" 

-ij 2 (v 2 ) t dxdt + (A ) / IVo^I 2 ^ 2 dxdt 

r 2 2 J R 



< 2 v 2 \V ip\ 2 dxdt 
Jr 
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Notice that \{ip 2 u 2 ) t = ipiptu 2 + uu t ip 2 . And use this fact by adding ip^u 2 to both 
sides: 

(2.24) / hip 2 u 2 ) t dxdt+{X-—) [ \V u\ 2 ip 2 dxdt < 2 / u 2 {\V ip\ 2 +ipipt) dxdt 
Jr 2 2 J R J R 

Now integrate this over R a = {B(x,r) xO<t<a<r} and drop the second 
integrand on the LHS. It is > 0. We'll choose a so that 

(2.25) / u 2 \t =(J dx > -max / u 2 dx. 

Jb(x,p') 2te(o,<T) Jb{x,t') 

This gives us from ( 12. 24ft : 

' ' 2 iP 2 ) t dxdt<2 [ ~ 2 "" ^ 2 

Ra * JRo 



(2.26) / -(wY),M<2 / n 2 (|V ^r + dxdt. 



Now consider the 



(2.27) max / u 2 dx < 2 / u 2 \ t=tT dx 

*e(0,cr) JB(x,r>) JB(x,r>) 

which we constructed above in (I2.25p . ip — 1 on R' and t is fixed, so 



max / u 2 dx < 2 u 2 \ t=a dx 

te(0,a) JB(x,r') J B(x,r>) 



< 2 u 2 i)\ =a dx 

JB(x,r') 



< 2 / u 2 ^ 2 \ t=a dx 

JB(x,r) 

Now recall (I2.26P integrated over R a . We will combine that with the | factor as 
well as the fact that the above is twice the (12.261) inequality. Putting these together 
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gives: 

(2.28) max / u 2 dx < 8A // u 2 (\ V ^| 2 + |V>V>t|) dxdt. 

te(0,<j) J B (x,r') JJBa 

Then use the non-negativity of the RHS and take the integral over R for 

(2.29) ~ 8 // m2 (I V o^! 2 + \Wt\)dxdt. 

Now we use a result on the bounds of the horizontal gradient of test functions in 
order to attain a constant independent of ip. Define ip(x,t) := ipi(t)ip2(^) i n the 
following way, using the cut-off functions of [10J : 



1 if t' < t < r 
if t < 



"02 0*0 



1 if d cc (x, 0) < r' 
if d cc (x, 0) > r 

and as before we have a bound on the horizontal gradient. This allows us to say 



r 2 1 

|2 i I „/.„/. I / / ^HG 1 



(2.30) |v ^r + \m < (r^% + — 

(r — r'Y r — v 

Which gives us 



te(o,°)Jnf '{r-r') 2 t-t"JJ r 



(2.31) max / u 2 dx < 8( ^ H< ^ + - — -) // v 2 c/xdt. 



This proves the result for (12.91) Now we return to ( 12.24ft 

/ i(V>V) t cfedt + (A - — ) f \V u\ 2 ^ 2 dxdt < 2 I m 2 (|V ^| 2 + ipipt) dxdt 
Jr2 2 J R J R 

and drop the first integrand. Integrated over B! gives 

(2.32) (A - [ |V «| V 2 dxdt <2 [ m 2 (|V ^| 2 + <M) dxdt 
2 Jr' Jr' 
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Again use the fact that ip = 1 on R' 

(A-—) / \V u\ 2 dxdt < 2 f u 2 {\V ip\ 2 + \ipip t \)dxdt 
2 Jr> Jri 

Which gives us 

C 2 f 

V u\ 2 dxdt < 2(A - — ) _1 / u 2 (| V ^| 2 + |<M|) dxdt. 



R> 2 j RI 

And using the previous bounds on ip, ip t and |Vo"0| 

(2.33) f \V u\ 2 dxdt < 2( \ -— + [ u 2 dxdt. 

Jr> 2 (r-r') 2 t-t> J R 

And the integral on the RHS is taken over R since it is non-negative. This proves 
the result in (|2.8|) . and so the corollary is proved. For the case p < p$ < we let 
v = u p l 2 be a weak supersolution to f)2.10p in R. We then have for <p = vif) 2 the 
following: 



and 



= X^ 2 ) = Xj(u p/2 ip 2 ) = -vf^XjUip 2 + 2^X j ^u p/2 
And we substitute these into ( 12.101) : 



m „ 

(v^ 2 )^/ 2 " 1 ^ dxdt+ V / A j (x,V u)(-u p / 2 ~ 1 X^ 2 + 2V'X j Vu p/2 ) rfxdt >0. 
2 . ^ 2 



which gives us 



/ (^ 2 )-w p/2 ~V <M + V / Aj(a;, V w)-w p/a ~ 1 Xj«V' 2 dxdt 
JR 2 J R 2 

>-2V / AjfaVoutyXjijM*/ 2 dxdt. 

7 = 1 ^ 
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At this point, notice that p < and the sign on the 2 in the RHS. We multiply by 
— 1 and proceed. 

J (^ 2 )M m p/2-i m< dxdt + J2 j Ajix^ou^uP^Xjiiij 2 dxdt 



m „ 

<2V / Aj(x, Vou)ipXjijju p/2 dxdt. 
j=i JR 

We calculate \{v t ) 2 = \u v l 2 u v l 2 ~ 1 u t and substitute into the first integrand: 

/ -^ 2 {v t ) 2 dxdt + jr I A j {x,V u)^u p/2 ~ 1 X J u^ 2 dxdt 
Jr 2 ~[ Jr 2 

<2V / Aj(x, W u)^X^u pl2 dxdt. 

Again we use the Cauchy and Young's inequalities on the RHS and recognize that 
u p = v 2 and V v = ^vP^XjU. 

f 1 m 

/ ~^ 2 (v t ) 2 dxdt + y~] e R tAj(x, V w)V W> 2 dxdt 

Jr 2 j=1 

<l I \Aj(x, V w)| 2 dxdt + 2 / v 2 \V^\ 2 dxdt. 
2 Jr Jr 

Now use the (S) bounds on A,- 

[ \i> 2 {v 2 ) t dxdt+\ [ \V v\ 2 ^ 2 dxdt < %- f |V u|V 2 dxdt+2 [ v 2 \V tp\ 2 dxdt 
Jr2 Jr 2 J R J R 

And from this point, the proof proceeds exactly like above. The case for < p < 1 
is modified in the following way: we must reverse t and use v := u p (x, —t). Then 
our test function must be = u p vp 2 . The proof proceeds in the same manner except 
there is a factor of x ^ 2p ~ 1 ' ) in front of the |Vof| 2 term. 



2.4 Nonlinear Iteration and bridge step 
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In order to run a Moser type iteration scheme, we needed the local bounds on 



(2.34) / \V u\ 2 dxdt <d f u 2 dxdt 

JR' Jr 



and 



(2.35) max / u 2 dxdt < C 2 / u 2 dxdt 

teloM J B (x,r') Jr 



Saloff-Coste shows in [32J that the parabolic Harnack inequality is equivalent to 
the doubling property 

fiB{x, 2R) < AfiB{x, R) 

and the Poincare inequality 



I \v-V B \ 2 dx < Cr 2 [ \Vv\ 2 dx 

J B(x,r) JB(x,r) 



' B(x,r) J B(x,r) 

as we have used in the proof above. This is done in the Riemannian manifold setting. 
He assumes the matrix of coefficients, ttjj ; £1X6 bounded and measurable, symmetric 
and elliptic. Saloff-Coste addresses the Moser iteration and on its dependence on 
both the Sobolev inequality and the Poincare inequality. Also, Grigory'an |15| 
showed the same for non-compact Riemannian manifolds. Now we will use the 
local bounds above to run the Moser iteration. 

Theorem 2.5. Let u be a positive, weak subsolution to 

m 
3=1 

in R with R' defined as above. If the structural hypotheses (S) and the general 
assumptions from the local bounds theorem hold, then for p > 1, there exists a 
constant 7 depending on X, C , n, and R such that 

(2.36) max «< 7(^1 1 1 u p dxdt) 1/p . 

Rl \R\ J Jr 
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The proof follows just as in the parabolic case. The constant 7 is slightly 
different. This is due to the fact that in the original parabolic case, the a,ij(x,t) 
are bounded and elliptic with constants A and h whereas the function Aj(x,^7ou) 
in the nonlinear case is bounded by a constant C and another A which is not an 
eigenvalue. Here 7 is the finite sum of terms 

(2.37) 7 = [O a (I)«8(^j + ^p) + m^ffy + 7^77) C B 

Proof: In order to prove the theorem, we use the lemma modeled after Moser 
Now for v G S 1,2 (R), define as before: 



(2 38) H ' Av):= \B^M\JJ x vdxdt 



(2.39) M r Jv) := — — max / v 2 dx 

\B r \ ie[0,r] Jn 



(2 - 40) ^'^^(o^iyj^ 1 ^ 

Lemma 2.6. Lei v G S^i?). Lei k = 1 + ^ for n G N, a = /3 = § T/ien, we 
/iawe /or some constants C\, C2 > 0: 

(2.41) H r ,, T (v k ) < d M r>T (v)n [C 2 H r ,, T (v) + D r ^ T (v 



which we proved above. In terms of H r>T (v), D rT (v) and M rr (v) we calculate 
our local bounds ( 12. 8 p and ( 12. 9 p as: 

(2.42) M rT (v) = max / v 2 dx 

\B r >\te(o,T) J B 



(2 43) - 8i ^- + ^mII/ dxdt 



( 2 - 44 ) ^ 8 ( j^Jf + t=F )h ' m 



(2-45) <8(^^ + -L 7 )^ r>T ( v ) 

Then for some C 3 depending on A, i? and n we have 

(2.46) M r ,»<C 3 # r ,>) 
And for 



D ' Av)= \B r x\ ,T')\ H™ 2dxdt 



we calculate: 



5 r x (0,r')| |(0, r)| v 2 y v ( r -r') 2 r - r' 

< 1 |(0 ' r)| 2fA C2 rv i 1 ) rr 

- |B r x(0,r)| 1(0, 7^)1 1 2 J \r-rT t-t' } JJ R 



v 

R 



< l 2 (A-^)-'(-^Lj + ^)// r >) 
r' 2 (r — r) z t — t 

gives for some C4 dependent on A, n and R 

(2.47) %'M<C 4 H r , T W 

So substituting these into the results of the lemma above, we get 

K, y {v k ) < C l CzH r , T {v) 2 -[C,H r , T (v) + C 2 H r , T (v)} 



= C 1 C 3 [C 4 + C 2 \H r , T (v)%H 1 , T (v) 



KG 



= 1 H r , T {v)~ +1 

(2.48) =lH r , T {v) k 
where 

7 = 4 CsCB (I)« [ (A-f)-^(^ + ^) + ^](S| B ,|(^ +7 2_))S 
We can repeat this exponentiation by k as in the parabolic case. Take 

(2.49) v v := W 2 

a positive subsolution in R where p v = pok u , p v > po > p' > 1 and k = 1 + 2/n, 
n E N is as above. 

Now we set up the domains for each of the iterations. Each iteration is defined 
in its own rectangle. 
Let 

(2.50) t„:=t 1 + T '" 



1 + TV 



and 



(2.51) r, = r l -^l 

1 + rv 



Let 

(2.52) K = B rv x (0,r,) 

Set 



(2.53) H u :=H rvTv {v v ) = T ± n ft v^dxdt. 

\ K v\ J J R v 



Now we will show the iteration. We notice that 
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(2.54) H v+1 = H rv+1Tv+1 (v v+1 ) = — ^ / f v^dxdt < ^H k u 

where 7 depends on R, n, A only. 

Now we will use this fact to achieve a maximum. First iterate down to Hq. Take 
any H u and notice that 

H v < 7 i^_! < 7^-1 



for some 72 > 7f . So 



H v < 7 2 ^-i 



Now use (H u _i) k to iterate again. 

= 72 (^-i)" < rtil^HUf 



— 72 n i/-2 



< 7 H + (^i) fe(7 (,-2)^ 3)fe2 



(i/) + (i/-l)fc+(i/-2)fc 2 rrk 3 
12 n v-Z 



. H + ^-l)fc+(^-2)fe 2 + ... + (l)fe"- 1 rrk" 
— 72 n 



< 73 ^0 



for some 



73 > 72 
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So we have 



which says that 



\R. 

Then taking k v roots: 



1 ff v PokV2 dxdt < 7 k^_L [I v Po/* dxdt ) 
M J Jr v \R\ J Jr 



(2.55) (t4-t// \v Po/ Ydxdt) l/k " ff v^ 2 dxdt. 

J Jr v \R\ J Jr 



letting v — )■ oo we conclude: 



maxM < 7 
i?' 



{ w\SL Mt) 



This proves Theorem 12.51 

Similarly we are able to achieve the inequalities 

(2.56) 1 ^\R\SS u ~ Pdxdt ) l/ P - m R [ F u 
and the crucial bridge step 

(2.57) W~\SS vP ° kV/2(ixdt - / / v- pokV,2 dxdt 

which are now seen to be independent of the nonlinear function Aj(x, Vqu). So the 
iterations for each of these inequalities may be run in a similar way with the slight 
difference that the constants will be altered in the manner explained above. 

2.5 Nonlinear Estimates on log(u) and Holder continuity 

Next we wish to obtain estimates on the function v := log{u) for a weak super- 
solutions u to 

III 

(2.58) u t = -J2 X jM x ^ u). 

3=1 
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If we are able to obtain estimates of the form 

(2.59) —J— / y/(v-C(x,r)) dxdt < N{v) 

\R 



R+ 



and 

(2.60) — }— / y/(C(x,r)-v) dxdt < N(v) 

\R I Jr- 

for some constant N depending on the function v, some constant C and rectangles 
R + and R~ subsets of R((x,t),r). Then we will be in a position to apply Aimar's 
[lj result for functions satisfying a BMO property with a lag mapping used above 
and complete the proof of the Harnack inequality. 

To begin, we will first let v = —log(u), = ip 2 ^ 1 , with ip > and of compact 
support on B(x,2r). We will let 



n 
u=l 



and 



1 if \x v — x°| < r 



if \x v — x v \ > 2r 

following [10J. Here we use (x?, x®, ...x° n ) to be the center of the ball in space. Then 
we calculate 

v t = {-log{u)) t = — 
u 

Xj<j> = Xji^u- 1 ) = 2i)Xji)u~ x - i?u- 2 XjU 

Xj(v) = Xj(-log(u)) = -Xjuiu- 1 ) 
And with (12.581) we have 
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ip 2 u~ l u t dxdt + 2^2 Aj^Xjtpu' 1 dxdt -J*] Aj^ 2 u~ 2 XjU dxdt > 0. 

j—l JJR j = i JJR 

Moving the second integrand to the RHS 

J J i) 2 u- x u t dxdt ~YljJ A ji )2u ' 2x 3 u dxdt > _2 J J A j^ x i^ u ~ l dxdt 
which gives us by substituting v t 

J J ^ 2Vt dXdt ^JJ A i^ 2u ~ 2X i u dxdt -~ 2 YljJ A 3^ X i^ u ~ l dxdt 
and hence 

J J i) 2 v t dxdt + Y2jJ Ajij} 2 u~ 2 XjU dxdt < 2 J J A^X^vT 1 dxdt 
Again we use the Cauchy and Young inequalities on the RHS to obtain 

J J ip 2 v t dxdt + J J Ajip 2 u~ 2 XjU dxdt 
- \lts II \Aj\ 2 u~ 2 ^j 2 dxdt + 2 II |V ^| 2 dxdt 



1 — JJR JJR 

Now use the boundedness of Aj(x, Vqu) from the (S) conditions on the second 
integrand on the LHS and on the first integrand of the RHS 

tp 2 v t dxdt + A / / tp 2 \ \7qu\ 2 u~ 2 dxdt 

R JJR 
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A 

< 
- 2 



J J |V w| 2 w"V 2 dxdt + 2 J J |V<y0| 2 dxdt 

At this point, notice that we have Xj(v) = —Xjuiu^ 1 ) from above which gives us 
|V u| = \Vqu\u~ 1 . 



J J ^ 2 v t dxdt + \ J J V 2 |V v| 2 dxdt < - J J \V Q v\ V dxdt + 2 J J \V ^\ 2 dxdt 
Gathering like terms 

(2.61) J I i?v t dxdt+(X --) J J 4) 2 \V v\ 2 dxdt < 2 J J |V ^| 2 dxdt 
Integrate from t\ to t 2 and observe that ip — 1 on B(x,r) is independent of t. 

[ v dx\\\ + (A - f * / |V(H 2 da;d* < 2(i 2 - h) [ |V Vf 

JB(x,r) 2 y tl J B {x,r) JB(x,r) 

Now we recall the Poincare lemma and the definition of Vb(£) = pj| / 



[Vb(t 2 )-Vb(ti)]|S(x,r)|+C7pr 2 (A-^) / ' / (« - Vi 

2 Jti JB(x,r) 

< 2(* 2 -fi) / |V ^| 2 cfe 

JB(x,r) 



,3 1 2 dx<it 



Divide by (t 2 — ^i) 



[^(t 2 ) - v B ( tl) ] + Cpr2(A _ a i f- r {v _ a dxdf 



(*2 - *l) ' 2 / (i 2 -t 1 ) 

< 2 / |V O -0| 2 dx 

JB{x,r) 

Take i 2 ->■ h 



ti JB(x,r) 
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dV B\r>f„ -M , n Jl(\ A \ I f.. t r \2 / o / ..|2 



, ,B{x,r)\+C P r\\- -) / (v - V B y dxdt < 2 / |V ^r dsc. 

2 Jb(x,t) J B{x,r) 

Divide through by \B(x,r)\ and use the bounds on |V(y0| < V r 



ll]B +t^-tt(\- 1 4) I (v-V B ) 2 dr,li < " ., , ;/?(.r./-i 



dt |S(x,r)| 2 7^^) J \B(x,r)\r 2 

We note that this inequality is the same form as (jl.l07p and that the proof is 
identical for both positive and negative t. This shows that v = —log(u) meets the 
BMO condition from Aimar [lj. So we may apply his results to finish the bridge 
step. Thus we have proved the Harnack inequality for positive, weak solutions to 

m 

ut = -^X*Aj{x, V w). 

3=1 

With the Harnack inequality in place, the Holder continuity with respect to the CC 
metric follows by the same argument as before following the pattern of J. Moser 
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